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CSeneral  Description 


Semiconductors  have  become  the  foxmdation  of  high  technology.  This  is  evidenced 
by  the  fact  that  every  piece  of  sophisticated  electronic  equipment  currently  contains  inte¬ 
grated  circuits  with  thousands  to  nearly  1.5  million  transistors.  The  improved  production 
of  such  equipment  is  thus  tied  up  with  the  effident  and  reliable  manufacture  of  very  small 
semiconductor  devices  and  with  accelerating  the  dedgn  of  mqga-transistor  drcuits.  These 
.  latter  issues  are  increasingly  tied  up  with  improved  modeling  and  better  understanding 
of  various  models  by  means  of  mathematical  and  computational  analysis.  This  is  because 
the  cost  of  experiments  is  very  high  and  because  many  of  the  phenomena  of  interest  are 
not  susceptible  to  direct  measurements.  The  subject  offers  a  uniqudy  rich  source  of  math¬ 
ematical  problems  such  as  free  boundary  problems  arising  from  modeling  surface  oxide 
growth  during  processing,  uniqueness  and  bifurcation  questions  for  systems  of  nonlinear 
elliptic  equations  in  the  case  of  the  drift-diffusion  model  of  a  angle  device,  existence  of  so¬ 
lutions  for  the  Boltzmaim  equation  in  more  realistic  device  models,  analysis  and  numerical 
methods  for  large  systems  of  differential-algebraic  systems  arising  in  circuitry,  and  so  on. 
New  challenges  in  modeling,  computations  and  analysis  have  recently  surfaced  due  to  the 
rapidly  decreasing  dimensions  of  devices. 


The  goal  of  this  program  was  to  foster  interaction  in  this  interdisciplinary  field  which 
involves  electrical  engineers,  computer  scientists,  semiconductor  physicists  end  mathemati- 
ciains,  from  both  university  and  industry.  The  program  particularly  encouraged  the  partic¬ 
ipation  of  numerical  and  mathematical  zmaJysts  with  backgrounds  in  ordinary  and  partial 
differential  equations  and  helped  get  them  involved  in  the  mathematical  aspects  of  semicon¬ 
ductor  models  and  circuits.  Leading  engineers  in  semiconductors  were  invited  to  present 
the  significeint  industrial  issues  as  well  as  to  concentrate  on  those  models  whidi  were  most 
relevant  to  mathematicians.  The  main  topics  of  the  program  were: 

I  or 


1.  Processing  modeling  (1  week)  July  15-19,  1991 

2.  Device  modehng  (l|  weeks)  July  22-31,  1991 


3.  Quantum  effects  (|  week)  July  Sl-August  2,  1991 

4.  Circuit  analysis  (1  week)  Aug.  5-9, 1991 

The  proceedings  will  be  published  as  two  volumes  in  the  series  IMA  Vol¬ 
umes  In  Mathematics  and  its  Applications  (Springer- Verlag).  • 

Week  1  •  Processing  Modeling,  July  15-19, 1991 

Pxocew  simulation  is  concerned  with  the  solution  of  physical  equations  which  model 
impurity  distribution,  oxide  growth  and  other  standard  wafer-processing  steps  (i.e.,  those 
necessary  to  fabricate  complete  device  structures)  such  as  thermal  annealing,  oxidation  and 
etching.  The  multitude  of  processing  steps  involves  a  corresponding  array  of  mathematical 
problems.  For  example,  stiff  nonlinear  sj^tem  of  parabolic  reaction-diffusion  equations  are 
used  to  model  point  defect  assisted  diffuskm;  oxidation  leads  to  singular  perturbation  pr  'b- 
lems  and  to  the  study  of  nonlinear  viscoelastic  time-dependent  equations  and  hyperbolic 
problems,  replete  with  shocks,  appear  in  the  study  of  reactive  ion  etchings.  Moreover,  in 
addition  to  the  traditional  numerical  analysis  questions  associated  with  the  above,  delicate 
computational  geometry  questions  arise  in  the  attempt  to  create  appropriate  moving  grids 
for  obtaining  the  impurity  distributions  within  domains  with  free  boundaries. 

Week  2  -  Device  Modeling  (l|  weeks),  July  22-31,  1991 

Device  modeling,  primarily  concerned  with  electrical  characterization  of  a  single  device, 
is  described  by  coupling  electromagnetics  -  Maxwdl  equations  or  some  of  thrir  approxi¬ 
mations  -  to  charge  carrier  transport  equations.  In  current  and  near  future  technolo^es, 
carriers  may  be  regar'^cd  as  pseudo  particles  and  a  generalized  Boltzmann  equation,  with  a 
momentiun-dependent  mass  and  a  complicated  collision  term  incorporating  a  fair  amount 
of  physics,  govern  the  transport  of  carriers.  Numerical  methods  which  are  based  on  the 
Monte  Carlo  approach  or  some  of  its  versions,  are  appropriate  for  understanding  the  basic 
transport  phenomena.  Simpler  models,  involving  the  first  few  moments  of  the  Boltzmann 
equation,  provide  the  tool  for  the  modeling  of  actual  devices,  with  all  of  their  complicated 
geometries  and  impurity  concentration  profiles.  The  traditional  model,  referred  to  as  the 
drift-diffusion  model,  comprises,  in  the  stationary  case,  a  system  of  three  nonlinear  equa¬ 
tions  of  elliptic  type,  and  has  received  most  of  the  attention  so  far  in  the  mathematical 
literature.  Although  existence  theorems  and  numerical  and  asymptotic  methods  for  this 
model  have  been  an  active  reseeirch  area,  some  basic  questions  such  as  uniqueness  and 
the  shape  of  the  I-V  curve  remain  imsolved  (from  the  viewpoint  of  rigorous  mathematics). 
Also,  numerical  approaches  to  three  dimensional  devices  have  not  yet  been  satisfactorily 
developed.  More  sophisticated  models,  which  fit  in  between  drift-  diffusion  and  Boltzmann, 
are  necessary  for  the  near  future  submicron  devices.  For  c.xample,  one  such  (hydrodynamic) 
model  in  2D  comprises  five  nonlinear  equations  of  a  mixed  hyperbolic-elliptic  type  even 


in  the  stationary  case.  Most  of  the  roc'themati'^  as  well  as  the  computational  aspects  of 
such  models  are  quite  open. 

Week  3  -  Quantum  Effects  week),  July  31  -  Aug.  2,  1991 

As  the  dimensions  of  modem  devices  decrease  to  ultrasmall  levels,  quantum  effects 
begin  to  play  a  major  role.  Examples  are  the  quantization  of  states  in  confined  regions  •  thin 
channels,  quantum-wires  etc.  -  and  tunneling  phenomena  in  resonant  diodes.  This  involves 
a  coupling  of  electromagnetics  directly  to  quantum  mechanics,  i.e.  to  eigenfunctions  of  the 
Schroedinger  equation.  New  transport  models,, such  as  the  Wigner  psuedo-differential 
equation  and/or  what  may  be  loosely  called  quantum-Boltzmann  have  to  be  applied  for  an 
adequate  description  of  "quantum  devices**.  Path  integrals,  psuedo-differential  operators 
and  numerical  techniques  for  eigenvalues  and  eigenfunctions  are  expected  to  provide  the 
mathematical  tools  for  this  area. 

Week  4  -  Circuit  Analysis,  August  5-9, 1991 

The  high  cost  of  integrated  circuit  fabrication,  both  in  time  and  materials,  forces  de¬ 
signers  to  simulate  their  complete  designs  in  as  much  detail  as  possible.  Typically,  the 
approach  of  choice  is  to  construct  a  differential  algebraic  (DAE)  system  that  describes 
the  circuit,  and  solve  that  system  numerically  for  a  ^ven  set  of  inputs.  More  specifically, 
the  circuit  is  assumed  to  be  an  interconnection  of  lumped  elements  (e.g.  resistors,  capaci¬ 
tor,  inductors  and  transistors),  characterized  by  a  relation  between  the  terminal  voltages 
and  terminal  currents.  The  equation  system  that  describes  such  a  network  is  then  con¬ 
structed  ULng  the  Kirchoff  voltage  and  current  laws,  and  this  leads  to  an  index  1  system  of 
differential-algebrmc  equations  where  the  state  variables  are  not  always  readily  identified. 
In  general,  the  system  of  equations  generated  from  most  circuits  is  very  stiff,  even  if  the 
algebraic  equations  are  not  present.  Also,  the  system  is  quite  nonlinear,  and  this  makes 
finding  an  initial  solution  which  satisfies  the  initial  conditions  and  the  algebraic  constrmnts 
rather  difficult. 

The  circuit  simulation  problem  can  be  broken  into  two  cl£isses:  solving  large  initial 
value  problems,  with  as  many  as  10,000  unknowns,  efficiently;  and  computing  periodic  or 
quasi-periodic  solutions  for  smaller  problems,  with  10  to  100  unknowns.  For  solving  the 
large  initial  value  differentieil-algebraic  systems,  the  topics  of  interest  are  techniques  for  stiff 
systems;  techniques  for  solving  large  nonlinear  systems  where  the  Jacobian  can  be  much 
more  than  one  rank  deficient;  multirate  integration  methods;  and  iterative  techniques  for 
very  nonsymmetric  problems.  For  computing  periodic  or  quasi-periodic  solutions  to  these 
DAE  systems,  topics  of  current  interest  are  existence  and  uniqueness  of  solutions,  high 
order  methods  for  very  stiff  problems,  spectral  methods,  mixed  time-frequency  methods 
and  other  types  of  collocation  methods.; 
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Week  1  •  Processing  Modeling,  July  15-19, 1991 
(Julian  Cole  will  be  in  residence  and  in  charge  of  week  1) 

Process  tiTmilmtinn  b  concerned  with  the  solution  of  physical  equations  which  model  impurity  distributiM, 
oxide  growth  and  other  standard  waler-procesaing  st^>s  (i.e.,  those  necessary  to  fabricate  complete  device 
structures)  such  as  thermal  annealing,  oxidation  and  etching.  The  multitude  dT  processing  steps  involvm  a 
corresponding  array  o(  mathematical  ptoblenos.  For  example,  stiff  nonlinear  system  parabolic  reactiM- 
diffusion  equations  are  used  to  model  point  ddect  assisted  diffusion;  oxidation  leads  to  singular  perturbatioo 
problems  and  to  the  study  of  nonlinear  viscoelastic  time-dependent  equations  and  hyperbolic  probfems, 
replete  with  shocks,  appear  in  the  study  of  reactive  km  etchings.  Moreover,  in  addition  to  the  tradition^ 
numerical  analysis  questions  associated  with  the  above,  delicaie  computational  geometry  questions  arise  in 
the  attempt  to  create  appropriate  moving  grids  for  obtaining  the  impurity  distributions  within  domains  with 
free  boundaries. 
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Monday,  July  IS 

UnloM  otbsrwiae  atated,  the  talk*  today  are  in  Conference  Hall  EE/CS  3*180 


9:00  am 

Hegbtraition  and  coffee 

Reception  Room  EE/CS  3-176 

9:30  am 

Welcome  and  Orientatbn 

Conference  Hall  EE/CS  3-180 

9:40  am 

Leonard  Bonicld 

.  Modeling  the  growth  and  annealing  of 

Motorola  Corporation 

dialocation  loopa 

Accounting  for  the  nuclcation,  growth  and  annealing  kinetics  of  dislocaUon  loops  is  beconung 
increasb^y  important  for  understanding  the  diffusion  behavior  of  the  dopants  that  are  commonly  used 
m  the  silicon  semiconductor  mdustry.  It  has  been  hypothesised,  for  example,  that  the  absorptioa  and 
emission  ot  silicon  bterstitials  by  end  of  range  dislocation  loops  has  a  signUicaat  transimt  effect  on  the 
diflusbn  of  iomimplanted  boron  during  rapid  thermal  processing  and  during  low  temperature  processing. 
A  better  understanding  of  this  transient  dHlusion  effect  is  desirable  because  it  wodd  b  turn  permit  better 
prediction  and  contrd  of  proceeses  used  to  create  some  dt  the  botoii*doped  structures  that  are  used  b 
current  technologies.  Altiiough  a  detailed  rinoulation  of  the  relevant  phenomena  is  not  currently  practical, 
something  can  be  learned  from  an  analysis  of  simpler  cases.  This  paper  examines  some  of  the  sdutiom  isnd 
approrimatioas  that  ate  b  the  liter^ure  for  the  preanire  Arid  around  a  pure  prismatic  iiriocttiou  bop. 
A  boundary  value  problem  is  formulated  for  the  growth  or  shrinkage  of  a  loop  due  to  Uie  bterae'don  of 
t'<.tetstitia]e  and  vacancies  rrith  the  preaeure  field.  An  issue  of  primary  concern  is  be  boundary  ccadition 
at  ibe  core  ci  the  bop.  Nunaerical  solutions  of  the  boundary  value  problem  are  presented  that  explore  the 
characteristics  of  the  model  and  their  relationship  to  the  observed  diffusion  phenomena. 

10:40  am  Coffee  Break  Reception  Room  EE/CS  3-176 

UKX)  am  P.  Lloyd  IC  technology  CAD  overview 

AT&T  BeU  Labs 

Ahiinet  This  paper  will  provide  an  overview  of  IC  Technology  CAD  b  terms  of  process  modeling,  device 
modeUng  and  circuit  analysis.  These  are  the  major  topics  of  the  1901  IMA  Summer  Program  on  Semiconduc¬ 
tors.  The  talk  will  discuss  the  btegratica  of  the  various  rimubtion  tods  and  discuss  bdustrial  applications 
in  process  developroent  and  product  design. 

2:00  pm  Martm  Giles  Pout  defect  diffusion  modeling  and  transient 

University  of  Michigan  diffusion  effects 

AMnet:  Dopant  diffusion  in  silicon  has  been  studied  for  several  decades,  but  is  still  the  subject  of  active 
research  and  unexpected  experimental  results  contuue  to  be  obtained.  Current  models  consider  di^ant  dif¬ 
fusion  to  occur  throu(^  interaction  with  pout  defects  -  silicon  vacancies  and  self-interstitials.  A  background 
concentration  of  pout  defects  exists  in  silicon  at  thermal  equilibrium,  and  they  can  also  be  generated  by 
processes  such  as  oxidation  and  ion  implantation,  and  can  be  removed  at  the  surface  or  at  dislocation  or  pre¬ 
cipitation  sites  withu  the  wafer..  Pout  defects  also  exist  m  various  charge  states  controlled  by  the  local  fermi 
level.  Macroscopic  diffusivity  is  the  net  result  of  dopant-defect  pair  formation  and  diffusion  for  each  charge 
state.  Part  of  the  difficulty  u  modeling  this  system  is  the  wide  range  of  timescales  involved,  from  electronic 
equilibrium  in  nanoseconds  through  dopant-defect  pair  formation  and  defect  diffusion  in  milliseconds  and 
seconds  to  macroscopic  dopant  diffusion  in  minutes  and  hours. 

For  ion  implanted  dopant,  the  process  of  ion  implantation  generates  point  defects  with  concentration  several 
orders  of  magnitude  above  the  background  level..  During  subsequent  .mnealing,  these  may  recombine,  r2q)idly 
diffuse  away,  or  may  coalesce  to  form  extended  defects.  While  excess  defects  remain,  dopant  diffusion  will  be 
greatly  enhanced.  Experimentally,  diffusion  is  seen  to  be  enhanced  by  factors  of  more  than  1000  for  times  on 
the  order  of  20  minutes  at  800  C  and  5  seconds  at  1000  C.  For  higher  implantation  doses,  extended  defects 
are  formed  in  a  narrow  band  in  the  tail  of  the  implantation  distribution.  Since  these  are  strong  sinks  for 
interstitials,  defect  concentrations  and  dopant  diffusivity  have  a  strong  variation  in  space  as  well  as  time. 
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During  this  talk  we  will  review  the  development  of  the  models  currently  used  for  diffusion  and  ion  im¬ 
plantation,  and  the  experimental  results  which  provide  the  physical  justification  for  the  noodels.  Examples 
will  illustrate  the  capabilities  offered  and  diffic^ties  encounte^  which  mi^^t  be  overcome  with  inproved 
mathematical  methods. 

4:00  pm  Vincent  Hall  802  IMA  Tea  (and  more!) 

(The  IMA  Lounge) 


X^Msday,  July  16 

Unless  otherwise  statedit  the  talks  today  are  in  Confiurence  Hall  EE/GS  S>180 
9:30  am  J.R.  King  Asymptotic  analysis  for  models  for  Impurity- 

University  of  Nottingham  de^t  pair  diffusion 

Aitinek  The  coupled  diffudon  of  impurities  and  pmnt  defects  in  nlioon  may  be  modelled  by  aystems  of 
reaction-diffusion  equations.  These  aystems  usually  contain  several  widely  different  time  and  q>aca  aeaks, 
and  asymptotic  methods  provide  a  systematic  means  of  redudng  their  oonplexity.  Yh  disoM  some  of 
these  physically  significant  asymptotic  limits,  identifying  the  apptopri^  parameter  ranges  in  each  case. 
Probl^  in  both  one  and  higher  dimensions  ve  conddwed.  Some  of  tiie  mathematically  unusual  features 
6t  the  resulting  simplified  problems  are  highlii^ted. 

10:30  am  Coffee  Break  Reception  Room  EE/CS  3-178 

lldX)  am  Ken  Suto  Atomic  dUfusbn  in  GaAs  with  contrrdled 

Tohdm  University  deviation  from  stoidiioroetiy 

AMnet:  Althou^  several  models  for  atomic  diffusion  in  GaAs  have  been  presented,  they  have  not  given 
strong  attention  on  the  effect  of  the  arsenic  v^>or  pressure,  i.e.,  the  deviation  from  atmchiometry,  or  some 
of  them  ate  thought  to  be  unrealistic. 

On  the  other  hand,  we  have  shown  that  the  crystal  growth  fitom  sdution  or  melt  under  applied  vapor  pressure 
i.e.,  temperature  difference  method  under  contt<dled  vapor  pressure,  can  be  explained  hf  the  equality  of  the 
arsenic  chemical  potentials,  and  the  dominating  point  dfiects  are  arsenic  interstitial  atoms  and  arsenic 
vacancies,  but  not  gallium  vacancies  and  gallium  intmstitials. 

On  the  basis  of  this  theory,  we  will  present  the  models  for  ^mic  diffusion  of  impurities  and  point  defects, 
particularly  with  attention  on  the  interstitial  atoms.  They  can  well  explain  the  arsenic  vapor  pressure 
dependence,  and  the  comparison  with  known  experiments  gives  reasonable  values  for  formation  energies  and 
migration  energies. 

Joint  work  with  Jun-ichi  Nishixawa. 

2:00  pm  Walter  Richardson  Reaction-diffusion  systems  arising  from 

University  of  Texas  at  San  AnUmio  phosphorous  diffusion 

Abstract:  A  five-species  reaction-diffusicm  system  is  derived  to  model  the  nonequilibrium  diffusion  of  phos¬ 
phorous  in  silicon  during  a  high-concealation  prediposition.  Mathematically  interesting  because  it  is  partly 
dissipative,  i.e.,  substitutionid  phosphorous  is  cor..'dered  immobile,  the  system  is  chailengmg  numerically 
because  it  is  stiff  and  ultimately  needs  to  be  solved  in  3  dimensions.  The  nunoerical  method  for  this  has 
been  used  to  implement  the  model  in  3-D  using  LSODP  and  various  matrix-free  iterative  methods  from  the 
NSPCG  package..  Results  indicate  that  a  large  Krylov  subspace  dimension  is  required  on  difficult  problems 
aud  that  for  realistic  goals  preconditioning  is  more  important  than  the  choice  of  an  iterative  method. 


Wednesday,  July  17 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 

9:30  am  Michael  J.  Johnson  An  interface  method  for  semiconductor  proces" 

IBM,  Lexington  simulation 
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AMrad:  The  diffusion  of  muiliple  species  of  dopants  in  silicon  at  high  temperatures  is  modeled  by  a  nonlinear 
parabolic  system  of  PDEs  ca  a  two>dimension^  regbn  with  a  moving  boundary.  A  numerical  aolutbu  umng 
the  L-stable  TRiaDF2  tune  integration  method  and  a  *bax  method”  spatial  discretisation  is  described. 
Details  are  given  of  the  methods  used  to  specify  curves,  to  noanipulate  cur/es,  and  to  define  arbitrary  simply- 
connected  redone  by  their  boundary  curves .  Numerical  experiments  comparing  alternative  chdces  of  matrix 
representation,  timestep  sise  selection,  and  operator  linearisation  are  described  and  analysed. 

10:30  am  Coffee  Break  Reception  Room  EE/CS  3-176 

IIKK)  am  P.  Dean  Gerber  The  extended  proximity  effect  equation  in 

IBM  electron  beam  lithography 

AMrad:  In  electron  beam  lithography  the  applied  dectrons  scatter  away  from  their  targeted  points  to 
nearby  (prosdmate)  points.  Hiis  is  known  as  the  proximity  effect,  and  correction  of  this  effect  is  necessary 
for  practical  lithogrm>hic  ^>pUcations.  Hie  scattering  is  commonly  modeled  by  a  multiple  diffusion  equation, 
so  the  emrection  problem  is  formulated  as  an  invene  diffunon  problem.  The  computational  complexity  of  this 
problem  is  a  functioo  of  the  number  of  sh^ies  in  the  lithographic  pattern,  which  can  run  in  the  hundreds  of 
millions.  We  dkcuas  this  problem,  which  is  both  ill-posed  and  ill-CMditioned,  and  obtain  a  practical  solution 
by  means  oi  a  “geometric  regularisation*  technique. 

3:00  pm  B.  Lojek  Rapid  thermal  annealing:  Diffusion  and  stress 

Motorola  Corporation 

AMrad:  Rapid  thermal  annealing  (RTA)  of  implanted  layers  is  a  new  technology  which  has  bem  widely 
studied  to  r^lace  conventional  thermal  processing.  However,  experimental  results  reveal  that  the  .physical 
phenomena  accompanybg  RTA  lead  to  many  anomalies  in  diffusion  behavior  which  in  the  literature  are 
usually  called  transient  diffusion.  Previously,  various  models  mostly  based  on  the  time  dependent  diffusivity 
or  epitaxial  realignment  of  damaged  layers  had  been  proposed  to  explain  these  phenomena  without  either 
solid  theoretical  support  or  unique  experimental  evidence.  There  is  no  theory  whidi  can  interpret  electrical 
and  chemical  profile  redistribution  during  RTA. 

It  has  been  found  that  the  prediction  of  the  lattice  position  and  the  electronic  configuration  of  the  implanted 
atoms  is  one  of  the  most  formidable  problems  m  diffusion  modeling.  Based  on  evaluation  of  available 
data  a  new  approach  it  presented.  The  impurity  atom  is  forced  by  the  local  tc^ological  constramt  of 
the  silicon  lattice  into  three  electronic  configurations:  for  group  III  inqmrities  -  acceptor  configuration 
with  tetrahedral  bonding  configuration  sp^;  natural  configuration  with  planar  sp*  configuration  and  finally 
the  nonbonding  (interstitial)  configuration  s’.  Similar  configurations  can  be  found  for  the  group  IV  of 
inqiurities:  donor,  natural  and  precipitates  configurations.  In  any  case  the  natural  electronic  configuration 
is  possible  only  if  a  vacancy  defect  is  in  the  vicinity  of  the  impurity.  More  important  however  is  that  the 
impurity  is  in  this  configuration  electrically  inactive.  For  example  in  this  context  only  the  VAs^  complex 
is  energetically  favorable  and  can  lead  to  electrical  deactivation.  Another  unwanted  side  effect  of  RTA 
processing  is  the  potential  to  generate  plastic  deformation  in  the  wafer.  During  plastic  deformation  a  large 
number  of  dislocations  is  introduced  into  certain  crystallographic  planes.  These  dislocations  have  impact 
not  only  on  impurity  diffusion  but  also  if  the  thermal  stress  is  large  enough,  permanent  crystallographic  slip 
occurs. 

The  main  purpose  of  this  paper  is  twofold:  first  to  review  the  RTA  processing  and  provide  evidence  that  the 
defect-impurity,  in  particular  vacancy-impurity  interaction  is  the  governing  process  which  determines  the 
electrical  activity  of  implant  and  subsequent  releixation  from  metastable  state  into  equilibrium.  Secondly, 
we  will  discuss  the  origin  of  thermoelastic  stress  originated  by  RTA  which  can  give  rise  to  the  slip. 


Thursday,  July  18 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 

9  30  am  Mark  Law  Automatic  grid  generation  for  diffusion  models 

University  of  Florida 


220 


Abitraci:  A  description  of  the  diffusion  equatbns  for  advanced  point-defect  based  models  of  dopant  diifusion 
in  silicon  will  be  presented.  These  modeb  require  the  solution  of  stiff,  non-linear  partial  differential  equations 
in  multiple  dimensions.  The  equations  require  accurate  solutions  over  time  ranges  from  milliseconds  to  hours, 
and  space  ranges  from  100  angstroms  to  500  microns.  In  addition,  these  problems  are  typically  solved  over 
domains  which  are  changing  in  time.  Consequently,  space  and  time  discretisation  is  critical  for  controlling 
both  errors  and  limiting  computer  time. 

Of  current  interest  is  manufacturing  simulation,  for  which  not  only  the  nominal  case  but  the  variances  ate 
solved  for.  For  this  situatbn,  control  the  spatial  and  temporal  discretisatbn  is  essential  so  that  the 
numerical  mor  can  be  controlled  to  be  less  than  the  process  variances  o(  mterest.  Techniques  for  saving 
these  problems  will  be  discussed. 


10:30  am 

Coffee  Break 

Reception  Room  EE/CS  3-176 

ll.-00am 

Andrsql  Strojwas 

Carnegie  Mellon  University 

New  numerical  techniques  for  efficient  simulaUon 
of  diffusion  and  lithography 

2K)0pm 

S.  Hamaguehi 

IBM,  Yorktown  Heights 

The  Boltsmann-Pokson  system  in  weakly 
eoUinonal  sheaths 

Abiirtet  The  energy  distribution  of  an  km  flux  bonobarding  a  cathode  surface  is  inoportant  information  u 
plasmaFpr  cessmg  took  for  etching  and  plasma-enhanced  dienaical  vapor  deposition  (PECVD).  In  particular, 
theoretical  understanding  of  such  bo  dynamics  in  recent  high-d^ity,  low-pressure  took,  in  which  plasrru 
sheaths  are  almost  collisionless,  k  of  significant  mterest.  In  thk  paper,  we  discuss  ion  kmetics  in  weddy 
oolkdonal  sheaths  described  by  boundary-value  probiems  of  the  Boltsmann-Poisson  system. 

First,  we  dkeuss  direct-current  (DC)  sheaths  or  steady-state  solutbns  of  the  system.  Assuming  that  elastb 
hard-sphere  collision  k  the  dominant  collkbn  mechanism  and  that  the  bn  mean-free  path  k  much  larger 
than  the  sheath  thickness,  we  obtain  the  ion  dktributbn  functbn  at  the  cathode  to  the  lowest  order  <d'  the 
coUisionality  parameter  (i.e.  the  ratio  of  the  sheath  thickness  to  the  bo  nrean  free  path).  Subsequently, 
angular  and  energy  distributions  of  the  bn  flux  are  calculated  based  on  the  dktribution  functbn,  whi^ 
provide  important  mformatbn  on  etching/PECVD  microscopic  profiles. 

Second,  we  discuM  radio-frequency  (RF)  sheaths,  where  a  boundary  value  of  the  system  (the  cathode  poten¬ 
tial)  k  a  periodic  function  in  time.  Since  the  RF  frequency  u  k  comparable  to  the  ion  plasma  frequency  utpi 
in  most  RF  processing  took,  the  oscillating  electric  field  exerts  a  nonlinear  force  (called  ponderomotive  force 
or  Milkr  force)  and  push  traveling  ions  mward  (toward  the  bulk  plasma)  b  an  RF  sheath.  Using  a  multi-time 
expansbn  method,  we  calculate  ponderoiimtive  force  m  the  limit  of  a  high  RF  frequency  ((w^j/w)’  <  1) 
and  estimate  decrease  of  the  bn  bombardment  energy  in  an  *RF  sheath  from  that  m  the  correspon^ng  DC 
sheath. 

The  theoretical  predictbns  described  above  are  also  compared  with  Monte-Carlo  simulations. 

5:30  pm  502  &  570  Vincent  Hall  Semiconductors  Pissa  Party 


Pizza,  salad  and  desert. 


Friday,  J  uiy  '» 9 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 
9:30  am  Contributed  Talks 

Contact  the  organizers  if  you  wish  to  give  a  contributed  talk  during  this  period. 

9'30  am  Sadasivan  Shankar  Numerical  solution  of  Boltzmann  equation  for 

University  of  Minnesota  electrons  in  glow  discharges  or  non-equilibrium 

plasmas 
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Absiraet:  With  the  advent  of  very  large  scale  mtegrated  circuits,  need  for  precise  control  of  various  micrO' 
electronic  devices  have  become  important.  The  study  of  plasma  is  an  analysis  of  the  dynamics  of  a  wide 
variety  of  particles  like  electrons,  positive  ions,  photons,  metastables,  free  radicals,  and  so  on.  The  talk  will 
be  divided  into  four  major  sections:  (i)  Background  of  plasma  processes  (ii)  Current  work  which  includes 
particle  approach  to  studying  plasmas,  (iii)  Contributions  and  conclusions  of  the  above  study,  and  (iv) 
Future  applicaticms  to  realistic  plasma  processes.  We  have  observed  that  in  nonequilibrium  plasmas  or  glow 
discharges,  due  to  low  degrees  ionisation  and  large  spatial  variations  in  electric  fields,  determination  of 
electron  behavior  is  dependent  cm  electron  velocity  distribution  function  (EVDF)  and  a  particle  iq>proach  is 
essential.  We  developed  an  approach  to  numerically  evaluate  the  EVDF  by  solving  the  l^ltsmann  equation 

electrons  in  a  plasma.  We  have  also  used  features  like  muHiiasking,  teetorixaiion,  and  opftmusfton  on 
supercomputers  like  Cray>XMP  and  Cray>2.  FVom  our  detaUed  study  of  electron  kinetics,  we  have  been 
able  to  calculate  transport  prc^erties  like  mesa  energy,  radial  and  longitudinal  fluxes,  and  denstfi'es.  Drijft 
ueloeiiies  were  also  computed  from  these  functions.  Tlie  output  of  the  algorithm  also  included  macroscopic 
transport  variables  like  diffusion  coefficients  and  mobilities.  These  properties  also  serve  to  validate  vvious 
modds  which  need  the  energy  dependence  of  transport  variables.  Evaluation  of  reaction  rate  constants  like 
ssctfsfisa  and  tonursfton  rife  ooasfsafs  are  necessary  to  understand  the  role  of  varbus  active  species  u 
plasma  reactors.  The  algorithm  supplies  these  as  funcUons  of  ebctron  energy  and  applied  electrical  fields. 
This  study  is  important,  sbce  it  is  useful  in  studying  the  chemistry  plasma  processes  and  the  link  between 
plasma  physics  and  plasma  processing. 

Joint  work  with  Klava  F.  Jensen. 

10:30  am  Coffee  Break  Reception  Room  EE/CS  3-176 


1  -  3:30  pm 


Industrial  Postdocs  Seminar 
Eleventh  Meeting 


The  format  of  the  seminar  is: 

1)  Presentatbn  of  projects  and  problems  from  industry  (3M  and  Honeywell)  on  which  the  bdustrial  postdocs 
are  working. 

2)  Informal  suggestions  and  discussion  among  the  participants. 

The  seminar  is  directed  by  Avner  Friedman  and  Walter  Littman.  The  topic  today  is  "Diffractive  Optics”. 
Those  wishing  to  participate  should  contact  A.  Friedman. 

The  SEMINAR  meets  in  Vincent  Hall  570 


Week  2  •  Device  Modeling  (1^  w«^ks),  July  22*31,  1991 
(FOiouk  Odeb  and  Bill  Coughran  will  be  in  residence  and  in  charge  of  weeks  2  and  3) 

Device  modelmg,  primarily  concerned  with  electrical  characterization  of  a  single  device,  is  described  by  cou¬ 
pling  electromagnetics  -  Maxwell  equations  or  some  of  their  approximations  -  to  charge  carrier  transport 
equations.  In  current  and  near  future  technologies,  carriers  may  be  regarded  as  pseudo  particles  and  a 
generalized  Boltzmann  equation,  with  a  momentum-dependent  mass  and  a  complicated  collision  term  in¬ 
corporating  a  fair  amount  of  physics,  govern  the  transport  of  carriers.  Numerical  methods  which  are  based 
on  the  Monte  Carlo  approach  or  some  of  its  versions,  are  appropriate  for  understanding  the  basic  transport 
phenomena.  Simpler  models,  involving  the  first  few  moments  of  the  Boltzmann  equation,  provide  the  tool  for 
the  modeling  of  actual  devices,  with  all  of  their  complicated  geometries  and  impurity  concentration  profiles. 
The  traditional  model,  referred  to  as  the  drift-diffusion  model,  comprises,  in  the  stationary  case,  a  system  of 
three  nonlinear  equations  of  elliptic  type,  and  has  received  most  of  the  attention  so  far  in  the  mathematical 
literature.-  Although  existence  theorems  and  numerical  and  asymptotic  methods  for  this  model  have  been  an 
active  research  area,  some  basic  questions  such  as  uniqueness  and  the  shape  of  the  1-V  curve  remain  unsolved 
(from  the  viewpoint  of  rigorous  mathematics).  Also,  numerical  approaches  to  three  dimensional  devices  have 
not  yet  been  satisfactorily  developed.  More  sophisticated  models,  which  fit  in  between  drift-  diffusion  and 
Boltzmann,  are  necessary  for  the  near  future  subniicron  devices.  For  example,  one  such  (hydrodynamic) 
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model  in  2D  comprises  five  nonlinear  equations  of  a  mixed  hyperbolic-elliptic  type  even  in  the  stationary 
case.  Most  of  the  mathematical  as  well  as  the  computational  aspects  of  8U<^  models  are  quite  open. 

‘  Monday,  July  22 

Cnloss  otherwise  stated,  the  talks  today  sure  in  Conference  Hall  EE/CS  3-180 

9:30  am  Giorgo  Baccarani  A  critical  review  of  the  fundamental 

University  di  B<dogna  semiconductor  equations 

Abrtnct:  Physical  models  for  numerical  device  simulation  are  reviewed  in  an  engineermg  perspective,  and 
their  derivation  from  more  general  physical  principles  is  examined.  More  specifically,  the  limitations  of  the 
basic  semiconductor  equations  within  both  the  drift-diffusion  and  the  hydrodynamic  models  are  discussed 
and  their  validity  range  is  assessed  with  emphasis  on  silicon  devices.  Special  attention  is  devoted  to  the  issue 
of  hot<arria  effects  and  related  sinmlation  problems.  Techniques  alternative  to  the  Monte  Carlo  ^>proach 
for  the  solution  of  the  Boltsmann  Transport  Equation  (BTE)  are  shortly  addressed.  Open  problems  both 
in  the  area  of  physical  models  and  numerical  te^niqnes  are  identified,  showing  the  need  of  fhrther  wo?k  to 
unprove  our  ability  of  reliably  predicting  device  perfcirmance  in  an  engineering  envir<»ment. 

Joint  work  with  Antonio  Gnudi  and  DaviJe  Ventura. 

10:30  am  Coffee  Break  Reception  Room  EE/CS  3-176 

2:00  pm  Peter  A.  Blakey  An  appUcat!onsK>riented  view  of  semiconductor 

Motorola  Corporaticm  device  simulation 

Abtinei:  This  presentation  will  focus  on  the  issues  faced  by  the  ievtlopen  and  asera  of  device  simulation 
codes.  An  undentanding  of  these  issues  can  help  specialised  mathematically-oriented  researchers  to  maximise 
the  utility  of  Uieir  efforts.  The  general  goals  of  the  talk  are:  to  supply  a  context  for  talks  that  deal  with 
specialised  topics;  to  provide  background  information  that  can  help  in  selection  and  prioritisation  of  research; 
and  to  give  nonengineers  an  interesting  glimpse  into  the  world  of  ^>plications. 

Code  developers  bridge  specialised  research  and  applications.  The  tasks  involved  in  implementing  a  device 
simulator  will  therefore  be  reviewed  first.  The  tasks  include:  selecting  physical  models;  selecting  numerical 
techniques;  selecting  a  software  engineering  strategy;  defining  capabilities  for  data  collection,  analysis  and 
visualisation;  and  coding.  Ttade-ofis  between  factors  such  as  generality,  efiiciency,  stability,  robustness,  and 
ease  of  implementation  are  required.  The  choices  in  each  area  will  be  surveyed.  A  by-product  of  this  exercise 
is  a  framework  for  classifying  device  simulators. 

The  second  part  of  the  presentation  will  focus  on  user  issues.  The  generic  modes  of  use  include:  obtaining 
physical  insight;  predicting  trends;  design;  optimisation;  developing  models  for  use  by  circuit  simulators;  and 
obtaining  parameters  required  by  circuit  models.  A  discussion  of  why  many  device  engineers  are  reluctant 
to  use  simulators  leads  to  the  identification  of  high  priority  research  areas.  Additional  discussion  will  center 
on  the  cost-benefit  concerns  of  management,  market  size  constraints,  and  manpower  limitations. 

The  last  part  of  the  talk  will  describe  specific  examples  of  the  successful  use  of  device  simulation  in  the 
development  of  microelectronic  products.  These  examples  will  include  microwave  diodes,  power  transistors, 
and  bipolar  and  MOS  technology  development. 

3:00  pm  Christian  Schmeiser  The  derivation  of  analytic  device  models  by 

TU-Wien-Austria  asymptotic  methods 

Abstract:  In  circuit  simulation,  device  models  should  be  as  simple  as  possible.  On  the  other  hand,  physic2dly 
sound  models  for  the  electrical  behaviour  of  semiconductor  devices  involve  nonlinear  systems  of  partial 
differential  equations  posed  on  domains  with  complicated  geometries.  Therefore  simplifications  have  to 
be  introduced  corresponding  to  certain  idealizing  assumptions.  By  the  use  of  asymptotic  methods  the 
simplification  procedure  can  be  carried  out  in  a  mathematically  justifiable  way. 

This  talk  gives  an  overview  of  recent  results  on  steady-state  voltage-current  characteristics  of  multidimen¬ 
sional  devices  as  well  as  a  new  approach  to  the  modelling  of  the  transient  behaviour  via  integral  equations 
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4:00  pm  Vincent  HaU  502 
(The  IMA  Lounge) 


IMA  Tea  (and  moK!) 


Tuesday,  July  23 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 

9:30  am  B.  Mmnershagen  On  the  consistent  modeling  of  impact  ionization 

University  of  Aachen  and  hot  electron  oxide  emission  in  NMOS 

devices 

Aistnict:  In  the  first  part  of  this  talk  experimental  results  on  hot  electron  oxide  emission  (gate  current) 
and  impact  ionization  (substrate  current)  are  reviewed  and  different  methods  of  modeling  these  effects  are 
studied.  This  leads  to  the  conclusion  that  neither  the  drift  diffusion  (DD)  model  nor  local  hot  electron 
models  (based  on  local  electric  field  and/or  electron  terr^erature)  can  be  used  for  a  consistent  rttodeling  of 
both  effects.  Moreover  it  is  shown  that  such  a  consistent  modeling  can  be  achieved  by  a  combination  of  a 
hydrodynamic  (HD)  model  and  a  nonlocal  field  line  based  lucky  electron  (LE)  model. 

In  the  second  part  of  the  talk  an  advanced  Monte  Carlo  (MC)  Model  is  used  u  a  physical  reference  to  confirm 
the  conclusions  drawn  at  first  from  experimental  results.  Finally  a  generalised  HD  model  that  is  consistent  to 
the  MC  model  is  introduced  and  it  is  demonstrated  that  the  internal  distributions  of  electrostatic  potential, 
electron  density,  electron  temperature  and  impact  ionisation  rate  resulting  firom  the  combination  of  the 
generalised  HD  model  and  the  nonlocal  LE  model  are  in  good  agreement  with  the  corresponding  results 
from  the  MC  model. 

10:30  am  Coffeo  Break  Reception  Room  EE/CS  3-176 

1  IKK)  am  Josef  Burgler  Combined  device-circuit  simulation  of  advanced 

ETH,  Zurich  semiconductor  structures 

Aisfrsef:  To  develop  and  optimise  semiconductor  devices  it  is  common  practice  to  advocate  numerical 
device  modeling  by  using  simulation  programs.  Typically  the  simulation  of  a  unit  is  split  into  levels,  starting 
with  process  simulation  (to  study  the  effect  ot  wafer-processing  steps)  and  device  simulation  (to  study  the 
electrical  behavior),  to  circuit  simulation  (to  study  the  electrical  features  of  an  ensemble  of  devices),  leading 
possibly  to  the  logic  simulation  of  a  whole  unit.  Unfortunately,  this  approach  nnakes  it  very  hard  to  model 
the  interaction  of  the  device  and  the  circuit  in  the  transient  case:  therefore,  in  order  to  optimise  the  transient 
behavior  of  a  particular  device  it  is  necessary  to  include  the  effects  of  both  the  drive  (i.e.,  control  part),  as 
well  as  the  load  circuit. 

During  the  last  year,  we  have  developed  a  software  environment  for  combined  device-circuit  simulation  studies 
dedicated  to  semiconductor  power  devices,  magnetic  field  sensws,  and  BiCMOS  structures  augmented  with 
simple  circuits.  In  our  work  we  discuss  some  critical  points  in  detail.  These  bclude  the  grid  generation  and 
adaptatbn,  the  scaling  of  the  unknowns  (i.e.  how  to  deal  with  high  voltages),  the  numerical  procedures  used 
to  solve  the  transient  problem  (with  special  attention  to  supercomputer  architectures),  the  assembly  of  the 
device  (element  assembly  versus  edge  assembly)  and  circuit  equations  as  well  as  the  physical  tnodeb  used. 

We  show  that  it  is  essential  to  chose  properly  the  numerical  techniques  to  obtain  accurate  and  reliable  results. 
Studies  of  the  turn-off  of  power  devices,  the  switching  of  BiCMOS  structures,  and  of  magnetic  field  sensors 
indicate  the  success  of  the  current  approach. 

2:00  pm  Carl  L.  Gardner  Electron  shock  waves  in  a  submicron 

Duke  University  semiconductor  device 

Abstract  The  hydrodynamic  model  treats  electron  flow  in  a  semiconductor  device  through  the  Euler  equa¬ 
tions  of  gas  dynamics  (with  heat  conuction)  coupled  to  Poisson’s  equation  for  the  electric  potential.  The 
hydrodynamic  model  PDEs  thus  have  hyperbolic,  parabolic,  auid  elliptic  modes. 

The  nonlinear  hyperbolic  modes  support  shock  waves.  Numerical  simulations  of  a  steady-state  electron 
shock  wave  in  a  submicron  seiruconductor  device  will  be  presented  using  a  steady-state  upwind  method.  The 
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electron  shock  wave  has  a  finite  width,  due  to  heat  conduction.  The  width  of  the  shock  profile  as  a  function 
of  heat  conduction  wili  be  analyzed. 


Wednesday,  July  24 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3*180 

9:30  am  Carlo  Cercignani  Methods  of  the  kinetic  theory  of  gases  relevant 

Politecnico  di  Milano  to  the  kinetic  models  for  semiconductors 

Ai$tnei:  When  the  transport  o^' charges  m  a  semiconductor  is  ccmsideted  (m  a  sufficiently  large  time  scale, 
then  the  motion  of  the  carriers  is  decidedly  influenced  hy  the  short  range  interactions  with  the  crystal 
lattice,  which  can  be  described,  in  a  classical  picture,  by  particle  collisions.  The  basic  tool,  in  this  situation, 
is  the  Boltsmann  equation,  which  may  be  simplified  by  excluding  the  short  range  interactions  between 
carriers,  which  only  play  a  role  when  the  particle  density  is  very  large.  The  basic  unknown  then  becomes 
the  distribution  function,  a  function  of  position  c,  velocity  v  and  time  t.  Even  if  we  exclude  the  interactions 
between  carriers,  we  obtain  a  nonlinear  equation  because  of  the  Pauli  exclusion  principle.  The  nonlinearity  is, 
however,  different  from  that  occurring  in  the  Boltzmann  equation  for  gases;  here  the  nonlwearity  can  make  life 
easier  because  it  ensures  that  /  is  bounded.  In  many  applications  the  nonlinear  term  can  be  neglected.  Then 
many  prop«ties  of  linear  transport  equations  can  be  applied.  There  ii,  of  course,  no  diflSculty  in  establishing 
existence  theorems  for  both  initial  and  boundary  value  problems.  As  an  approximate  computational  tool  it 
is  important  to  recall  the  variational  principle  for  the  ste^y  case  in  a  given  electric  field,  which  can  be  useful 
to  equate  i^obal  quantities  such  as  currents  and  jumps  near  contacts  and  interfaces  in  an  accurate  way. 
One  can  also  study  the  Green's  function  and  particular  solutions  that  are  extremely  useful  in  order  to  gain 
insight  into  the  qualitative  behavior  of  solutions.  The  talk  will  be  devoted  to  a  survey  of  the  abovementioned 
topics. 

10:30  am  Coffee  Break  Reception  Room  EE/CS  3-176 

11:00  am  F.  Poupaud  Boundary  value  problems  in  semiconductors 

Univeraitd  de  Nice  for  the  stationary  Vlaaov-Maxwell-Boltsmann 

equations 

Absinei:  The  works  of  R.J.  Diperna  and  P.L.  Lions  [1]  have  allowed  significant  progress  in  the  study 
of  Cauchy  problems  in  free  spaces  for  kinetic  models.  However,  there  are  few  results  on  boundary  value 
problems  in  this  field.  Here  we  present  some  recent  results  on  the  Vlasov  Maxwell  system  that  provide 
stationary  solutions  of  boundary  value  problems  for  some  kinetic  models  of  semiconductors.  The  study  of 
stationary  solutions  of  transport  equations  for  charged  particles  in  selfconsistent  fields  has  been  initiated  by 
C.  Greengard  and  P.A.  Raviart  [2].  They  have  given  an  analysis  of  the  Vlasov-Poisson  system  in  a  one- 
dimensicmal  geontetry.  The  techniques  developed  by  the  author  [3,4]  to  construct  stationary  solutions  for 
kinetic  modeb  in  any  kind  of  geometry  are  quiet  different  of  those  of  [1,2].  They  are  based  on  the  use  of 
upper  solutions  of  the  Vlasov  equation.  In  the  case  of  repulsive  interparticular  forces,  this  method  provides 
a-priori  estimates  that  allows  one  to  apply  the  Schauder  fixed  point  theorem.  Up  to  now  this  analysis  takes 
into  account  only  the  linear  collision  phenomena,  that  is  to  say  interactions  of  charged  particles  (electrons 
and  holes)  with  phonons  and  impurities  in  a  non  degenerate  semiconductor.  Some  works  are  in  progress  for 
the  analysis  of  non  linear  Boltzmann  opeiators  modeling  degeneracy  effects  and  binary  collisions. 

References 

1)  R.J.  D'nerna  and  P,L.  Lions,  Global  weak  solutions  of  Vlasov-Maxwell  systems,  Comm.  Pure  Appl. 
Math.,  XLII,  1989,  pp.  729-757 

2)  C.  Greengard  and  P.A.  Raviart,  A  boundary  value  problem  for  the  stationary  Vlasov-Poisson  system  : 
the  plane  diode.  Comm.  Pure  Appl.  Math.,  VLIII,  1990,  pp.  473-507 

3)  F.  Poupaud,  Solutions  stationnaires  des  equations  de  Vlasov-Poisson,  C.  R.  Acad.  Sci.  Paris  serie  1,  311,: 
1990,  pp.  307-312 

4)  F  Poupaud,  Boundary  value  problems  for  the  stationary  Vlasov-Mcixwell  system.  Forum  N'ath.,  to  appear 
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2:00  pm  Luis  G.  Reyna  On  the  treatment  of  the  collision  operator  for 

IBM  Watson  Research  Center  hydrodynwiic  models 

Abstract  The  study  of  the  charge  transport  in  semiconductor  devices  has  been  traditionally  based  on  the 
Drift-Diffusicm  model.  Smaller  devices  need  to  take  into  account  ballistic  electrons  effects.  An  almost  com¬ 
plete  description  of  these  effects  can  be  obtained  form  computationally  expensive  Monte  Carlo  simulations. 
An  alternative  a4;>proach  is  given  by  Hydrodynamic  models. 

In  this  work  we  propose  an  alternative  treatment  of  the  collision  operator  in  the  Boltcman  equation  leading 
to  Hydrodynamic  models.  We  start  with  a  trial  displaced  Maxwellian  function  for  the  distribution  function 
and  by  explicit  integration  of  the  collision  term,  we  obtain  a  set  of  differential  equations  that  lies  in  between 
the  Drift-IMffusion  model  and  the  Hydrodynamic  model. 

We  present  numerical  results  obtained  from  these  equations  for  a  simple  one  dimensional  N+  N  N+  structure. 
We  include  the  interaction  of  electrons  with  acoustic  uid  optical  phonons  and  discuss  the  poasibiEty  of 
including:  non-parabolic  bands  and  intervalley  scattering.  The  present  approach  could  be  suitable  for 
nmulatisns  baaed  on  a  mixed  Monte  Cark>-Hydrodynamic  description. 

Joint  work  with  Andres  Saul. 


Thursday,  July  25 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  HaU  Efi/OS  3*180 

9:30  am  Thomas  Kerkhoven  Numericdly  stable  three  dimensional  finite 

University  of  Illinois,  Urbans  element  discretisnAions  of  drift-diffusion 

equations 

Aksfrsct:  Hie  drift-diffusion  current  continuity  semiconductor  equations  ate  usually  brought  in  self-adjoint 
form  in  terms  of  Slotboom  variables  by  an  exponentiiJ  upwinding  technique.  The  equations  for  the  Slotboom 
variables  ve  usualy  diacretised  by  the  well-established  box  method.  This  box  method  can  be  considered  as 
a  Petrov-Galerkin  method  for  a  piecewise  linear  approximatimi  with  piecewise  constant  test-functions.  On 
the  Delaunay  triangulatbn  in  two  dimenrions  the  commonly  employed  box  method  discretisation  yields  a 
discretised  system  with  a  maximum  principle  for  the  Slotboom  variables. 

Because  the  piecewise  constant  test-functions  are  not  of  regularity  complications  arise  in  applying  stan¬ 
dard  finite  element  err<»  analysis  to  this  Petrov-Galerkin  method.  These  complications  ^  be  circumvented 
by  the  foUowmg  observations.  The  drift-diffusion  nnodel  for  the  Sbtboom  variables  is  described  by  an  elliptic 
gradbnt  equatbn..  Through  Gauss’s  law  we  show  that  the  Petrov-Galerkin  discretisatbn  of  an  elliptic  gra¬ 
dient  equatbn  for  a  pbcewise  Unear  approximation  function  depends  only  on  the  size  of  the  mtegrals  of  the 
test  functions  over  the  faces  of  the  elements,  and  on  the  regularity  of  the  coefficients  in  the  elUptic  gradient 
equation.  We  demonstrate  that  this  observation  impUes  immediately  that  m  two  dimensions  a  box  method 
in  which  the  boxes  bisect  the  finite  element  faces  is  equivalent  to  the  pbcewise  Unear  finite  element  method. 
Hence,  regularity  complications  can  be  avoided. 

However,  in  three  dimensions  a  box  method  is  only  equivalent  to  the  piecewise  Imear  finite  element  db- 
cretization  if  the  boxes  partition  the  faces  of  the  elements  b  three  equal  parts.  This  impUes  that  the  box 
method  based  on  perpendicular  bisecting  planes  and  the  piecewise  linear  finite  element  method  need  not 
generate  identical  stress  matrices  for  the  Laplacean  in  three  dimensions.. 

However,  the  box  method  based  on  perpendicular  bbectors  yields  a  maximum  principle  on  a  Delaunay 
tetrahedryzation  i/hereas  the  piecewise  Unear  finite  element  method  r  eed  not.  We  introduce  an  alternative 
Petrov-Galerkin  method  which  is  equivalent  to  the  box  method  based  on  perpendicular  bbecting  planes  but 
which  circumvents  problems  with  reduced  regularity.. 

10:30  am  Coffee  Break  Reception  Room  EE/CS  3-176 

1100  am  M.  Rudan  Integrated  toob  for  device  design 

University  of  Bologna 
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Absirxici:  Process  and  device  simulators  have  found  a  wide  acceptance  among  designers;  these  tools  are 
providing  reliable  results  that  allow  one  to  cut  the  cost  associated  with  the  laboratory  experiments  required 
to  investigate  different  tradeoffs.  Even  though  the  software  environment  to  support  device  design  has  had  a 
steady  improvement  in  the  last  few  years,  the  tuning  of  the  process  parameters  is  usually  left  to  the  process 
designer  who,  on  his  own,  has  to  choose  the  tradeoff  among  many  objectives  in  competition  with  one  another. 
From  this  point  of  view,  the  availability  of  an  integrated  synthesis  system  that  supports  an  optimized  design 
of  scmicouductor  devices  is  likely  to  be  an  useful  tool  for  a  timely  and  high-quality  design  development. 

The  first  part  of  the  talk  will  outline  the  development  of  an  optimization  tool  based  on  our  two-dimensional 
device-analysis  code  HFIELDS.  The  activity  is  carried  out  in  the  frame  of  an  EEC-sponsored  Project,  STORM. 
The  main  objective  of  the  optimization  part  of  the  Project  is  the  development  of  an  automatic  optimizer  for 
process  and  device  design,  whose  goal  is  computing  optimal  values  for  the  process  parameters  according  to 
given  targets  and  constraints. 

In  the  second  part,  a  brief  description  is  given  of  how  the  optical-generation  phenomena  have  been  incorpo¬ 
rated  into  HFIELDS.  This  goal  has  been  achieved  by  introducing  a  number  of  optical  windows  and  interleaved 
matx^rial  layers  through  which  a  radiation  with  arbitrary  spectrum,  incidence  angle,  and  polarization  state 
enters  the  crystal,  and  by  evaluating  the  corresponding  generation  rate  at  each  node  of  the  discretisation 
grid.  The  code  equipped  with  this  new  capability  makes  the  description  of  realistic  semiconductor  optical 
sensors  feasible. 

Joint  work  with  M.  C.  Veccbi,  G.  Verzellesi  and  Zs.  M.  Kovacs. 

2:00  pm  Henry  H.K.  Tang  Transport  theoretic  foundation  of  generalized 

IBM,  East  Fishkill  hydrodynamics 

Ahsirtci:  Some  of  the  central  issues  of  modeling  submicron  devices,  both  theoretical  and  computational,  are 
examined  from  the  viewpoint  of  two  complementary  approaches.  1.  A  macroscopic  theory  of  generalized 
hydrodynamics  is  derived  from  a  Boltzmann-Vlasov-type  (BV)  transport  equation.  2.  A  finite-difference 
method  is  devebped  to  solve  the  time-dependent  BV  equation. 

The  non-equilibrium  dynamics  of  hot  electrons  in  strong  fields  is  studied  from  the  framework  of  a  micro¬ 
scopic  BV  equation.  We  reformulate  the  cdlision  integrals  of  the  BV  equation  by  means  of  a  relaxation-time 
representation.  Each  type  o(  integral  -  due  to  electron-electron,  electron-phonon  or  electron-impurity  scat¬ 
tering  -  is  rewritten  as  a  series  of  relaxation- time-like  terms,  analogous  to  a  partial- wave  decomposition  of 
the  underlying  scattering  matrix.  IVuncation  of  these  series  leads  to  a  class  of  generalized  relaxation-time 
models  which  incorporate  much  of  the  salient  features  of  the  original  collision  integrals,  but  which  can  be 
more  readily  solved. 

Generalized  hydrodynamics  is  formulated  in  terms  of  the  hierarchy  of  moment  equations  of  the  BV  equation. 
The  simplest  version  consists  of  a  continuity  equation,  a  generalized  momentum  equation  and  equations 
of  the  temperature  tensor,  ^t  reduces  to  the  standard  drift-diffusion  model  and  classical  hydrodynamics 
under  the  restrictive  assumpUon  of  local  equilibrium  of  the  electrons.  We  emphasize  the  significance  of 
the  temperature  tensor  field  which  characterizes  non-equilibrium.  The  fundamental  problems  of  closing  the 
moment  equations  and  computing  the  transport  coefficients  are  transparent  and  can  be  solved  readily  in  a 
class  of  relaxation-time  models.  Our  approach  is  novel  in  several  respects.  First,  the  essential  physics  of 
microscopic  non-equilibrium  processes  is  included.  Second,  no  explicit  assumption  about  short  mean  free 
paths  of  the  electrons  is  invoked.  The  formalism  applies  even  in  the  ballistic  regimes  where  the  collisional 
frequency  is  low..  Third,  our  analysis  leads  to  new  computational  schemes  for  the  intrinsically  non-local 
transport  parameters  (such  as  mobility  and  thermal  conductivity),  which  go  beyond  the  confines  of  the 
classic  Hilbert-Chapman-Enskog-type  methods. 

Many  important  physics  i‘.sues  of  electron  transport  over  small  spatial  dimensions  and  time  scales  require 
detailed  information  of  the  electron  distribution,  A  finite-difference  technique  is  developed  to  solve  the  time- 
dependent  BV  equation.  This  method  is  based  on  an  accurate  and  efficient  solver  which  solves  the  collisionless 
Boltzmann- Vlasov  equation.  In  the  general  problem,  with  the  help  of  the  relaxation-time  representation, 
the  collision  integrals  are  expressed  as  simple  integral  operators  over  the  momentum  space.  Regarded  as 
'‘pseudo-potentials  ’ ,  the  collision  terms  (on  the  right  side  of  the  BV  equation)  and  the  self-consisteiit  coulomb 
field  (on  the  left  side)  are  treated  numerically  on  an  equal  footing.  These  ideas  are  tested  in  a  transport 
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model  in  I'dimensional  coordinate  space  and  3-diinensional  momentum  space.  We  discuss  the  prospects  of 
applying  this  scheme  to  problems  in  2-dimen8ional  coordinate  space  and  for  realistic  devices. 

?-3:30  pm  Emad  Fatemi  Simulation  of  n’*‘nn+  channel  using  finite 

IMA  difference  techniques  for  Boltzman  equation 

Abstract:  In  this  talk  we  discuss  simulation  of  a  channel  using  finite  difference  discretization  of 

the  Boltzman-Poisson  system  in  lD(space)*lD( velocity).  We  use  a  first  order  upwind  method.  The  collision 
operator  is  approximated  by  relaxation  to  a  Maxwellian.  The  simulation  is  compared  with  the  Hydrodynamic 
model.  Also  a  numerical  investigation  is  done  with  respect  to  the  heat  conduction  term  in  the  Hydro  model. 

5:30  pm  Semiconductors  BBQ  Picnic  Courtyard  between  Vincent  and  Murphy  Halls 


Fridc^,  July  26 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  HaU  EE/CS  3-180 
d:30  am  P.  Degond  On  the  Child-Langmuir  law  for  semiconductors 

ENS-Cachan,  France 

Abstract:  The  Child-Langmuir  injection  law  is  widely  used  in  vacuum  diodes  modellmg.  It  describes  the 
emission  of  an  electron  beam  at  a  cathode  when  the  energy  of  the  emitted  particles  is  small  compared  with 
the  applied  potential.  In  the  early  eighties,  the  physicists  M.S.  Shur  and  L.F.  Eastman  proposed  a  model 
of  biJlistic  transistor  based  on  the  Child-Langmuir  model  Such  models  can  be  mathematically  obtained 
by  means  of  an  asymptotic  analysis  of  a  singular  perturbation  problem  for  the  stationary  Vlasov-Poisson 
equation  (in  the  vaccuum  diode  case)  or  the  Vlasov-Poisson-Boltsmann  equation  (in  the  semiconductor  case). 
These  modek  produce  an  accurate  description  of  the  behaviour  of  the  carriers  close  to  the  cathode  contacts 
and  can  be  used  to  design  precise  boundary  conditions  especially  for  particle  or  Monte-Carlo  calculations. 

10:30  am  Coffee  Break  Reception  Hoorn  EE/CS  3-176 

11:00  am  Michael  Sever  Symmetric  forms  of  energy-momentum  transport 

Hebrew  University  models 

Aisfnicf’  The  energy-momentum  transport  models  of  carrier  flow  in  a  semiconductor  admit  entropy  functions 
and  thus  a  generalized  symmetric  form.  A  simpler  such  form  can  be  obtained  for  the  reduced  ^ass- 
momentum*^  system,  thinking  of  the  energy  equation  as  essentially  determining  the  carrier  temperature 
distribution.  The  use  of  such  a  symmetric  form  leads  naturally  to  suitable  choices  for  boundary  conditions, 
discretization  schemed,  and  regularization  when  discontinuous  solutions  appear. 

2:00  pm  Joseph  W.  Jerome  Energy  models  for  one-carrier  transport  in 

Northwestern  University  semiconductor  devices 

Abstract  Moment  models  of  carrier  transport,  derived  from  the  Boltzmann  equation,  have  made  possible 
the  simulation  of  certain  key  effects  through  such  realistic  assumptions  as  energy  dependent  mobility  func¬ 
tions.  This  type  of  global  dependence  permits  the  observation  of  velocity  overshoot  in  the  vicinity  of  device 
junctions,  not  discerned  via  classical  drift-diffusion  models,  which  are  primarily  local  in  nature.  It  has  been 
found  that  a  critical  role  is  played  in  the  hydrodynamic  model  by  the  heat  conduction  term.  When  ignored, 
the  overshoot  is  inappropriately  deimped.  When  the  standard  choice  of  the  Wiedemann-Franz  law  is  made 
for  the  conductivity,  spurious  overshoot  is  observed.  Agreement  with  Monte-Carlo  simulation  in  this  regime 
has  required  empirical  modification  of  this  law,  as  observed  by  IBM  researchers,  or  nonstandard  choices.  In 
this  lecture,  simulations  of  the  hydrodynamic  model  in  two  dimensions,  as  well  as  simulations  of  a  newly 
developed  energy  model,  the  RT  model,  will  be  presented.  These  simulations  have  been  carried  out  by  Chi- 
Wang  Shu  of  Brown  University.  The  RT  model  intermediate  between  the  hydrodynamic  and  drift-diffusion 
model,  was  developed  at  the  University  of  Illinois  to  eliminate  the  parabolic  energy  band  assumption,  and 
to  reduce  the  spurious  overshoot  with  physically  consistent  assumptions.  The  algorithms  employed  for  both 
models  are  the  essentially  nonoscillatory  shock  capturing  algorithms, >  developed  at  UCLA  during  the  last 
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decade.  Some  mathematical  results  will  be  presented,  and  contrasted  with  the  highly  aeveloped  state  of  the 
drift-diffusion  model. 


Monday^  July  29 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 


9:30  am 

Michael  Ward 

Stanford  University 

Some  examples  of  singular  perturbation  methods 
in  device  modeling 

10:30  am 

Coffee  Break 

Reception  Room  EE/CS  3-176 

11:00  am 

William  Coughran 

AT&T  BeU  Labs 

The  role  of  numerical  algorithms  and 
visualization  in  complex  semiconductor 
simulation 

2:00  pm 

Herbert  S.  Bennett 

NIST 

Physics  for  device  simulations  and  experimental 
verification 

The  motivations  f<^  using  computers  to  simulate  the  electrical  characteristics  of  transistors  are 
discussed.  Our  work  and  that  of  others  in  the  area  of  device  physics  and  modeling  are  described.  We 
compare  conventional  device  physics  with  an  alternative  approach  to  device  physics  that  is  more  directly 
traceable  to  quantum-mechanical  concepts.  We  then  apply  this  new  approach  to  quasi  neutral  regions, 
space-charge  regions,  and  regions  with  high  levels  of  carrier  injection.  Examples  of  applying  quantum- 
mechanically-based  device  ph^ics  to  energy  band  diagrams  for  bipolar  transistors  are  given.  The  limits  for 
using  theoretical  results  from  uniform  media  in  numerical  simulations  of  devices  with  large  concentraticm 
gradients  are  discussed.  Calculatbns  ot  the  effective  intrinsic  carrier  concentrations  for  gallium  arsenide  and 
silicon  are  also  given  along  with  published  data.  In  addition,  calculations  of  the  mobilities  for  GaAs  that 
are  baaed  in  part  on  quantum-mechanical  phase  shifts  are  compared  with  published  data.  We  then  conclude 
with  a  discussion  ct  the  requirements  for  verifying  and  calibrating  device  simulators  for  the  submicrometer 
domain. 

Joint  work  with  Jeremiah  R.  Lowney. 

3-3:30  pm  Tom  Seidman  Mathematical  treatment  of  the  time-dependent 

U.  of  Maryland,  Baltimore  Cty.  drift/diffusion  model 

Absimei:  The  drift/diffusion  mode!  is  considered  under  general  enough  conditions  to  include,  e.g.,  ^standard** 
impact  ionization  source  terms.  Existence  is  shown  for  the  general  model  and  then,  under  somewhat  more 
restrictive  conditions,  a  constructive  (contractive  iteration)  argument  gives  both  existence  and  uniqueness. 


Tuesday,  July  30 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 
10:30  am  Coffee  Break  Reception  Room  EE/CS  3-176 

9:30  am  Thomas  Kerkhoven  Numerical  simulation  of  quantum  wires  in 

University  of  Illinois,  Urbana  periodic  heterojunction  structures 

Abstract:  Electron  confinement  in  quantum  planes,  wires,  and  dots  is  modeled  by  Schrodinger’s  equation 
coupled  with  Poisson’s  equation  for  the  electrostatic  potential.  A  self-consistent  solution  to  this  system  of 
equations  is  obtained  as  follows.,  Schrodinger’s  equation  is  solved  by  subspace  iteration.  Poisson’s  equation 
is  solve  d  by  a  suitably  preconditioned  conjugate  gradient  algorithm.  Self  consistency  of  ihe  coupled  model 
is  obtained  by  a  Jacobian-free  implementation  of  Newton’s  method. 
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The  periodic  structure  can  only  be  employed  for  the  design  of  quantum  wires  if  the  following  conditbns  can 
be  met.  The  periodic  structure  allows  a  ground  state  wave  function  which  b  bound.  The  occupancy  of  the 
non-bound  states  can  be  made  neglipble  with  respect  to  the  occupancy  of  the  bound  ground  state.  This 
second  conditimi  requires  again  that  the  separation  between  the  bound  state  and  the  lowest  non  bound  state 
is  sufficient. 

We  briefly  describe  our  algorithm  and  proceed  to  discuss  in  detail  the  design  of  the  periodic  heterojunction 
structure  which  satisfies  the  requirements  mentioned  above. 

11:00  am  Peter  Szmolyan  Analysis  of  the  Gunn  effect 

Technical  University  Vienna 

Ah$iract:  If  a  constant  voltage  above  a  certain  threshold  is  applied  to  a  piece  of  semiconductor  maierial  with 
negative  differential  resistance  periodic  current  oscillations  are  observed.  The  current  peaks  are  due  to  dipole 
waves  which  are  generated  periodically  at  one  contact  be  device  and  leave  at  the  other  contact.  We  give  a 
refined  analysis  of  the  classical  explanation  of  the  Gunn  effect  as  traveling  waves  on  an  infinite  domain.  We 
show  that  under  appropriate  boundary  conditions  periodic  solutions  on  a  finite  domain  are  generated  by  a 
Hopf  bifiucation.  We  give  an  asymptotic  analysis  of  the  periodic  generation,  propagation,  and  destruction 
of  the  dipole  waves. 

2:00  pm  Patrick  S.  Hagan  High  field  semicmiductor  equations 

Los  Alamos  National  Lab 

Abstract:  Semiconductors  are  commmily  modeled  and  analysed  at  two  different  levels  of  description.  At  the 
macroscqpic  level,  the  Basic  Semiconductor  Equations  (BSEs)  are  a  system  of  diffusbn-drift  equations  which 
describe  the  densitites  of  electrons  and  holes  in  real  space.  At  a  much  more  fundamental  level,  Boltzmann 
equations  describe  the  densities  of  electrons  and  hdes  m  phase-space.  Although  the  BSEs  ate  quite  a^i'urate 
when  ebctric  forces  are  not  too  large,  it  is  well-known  that  the  BSEs  become  increasingly  inaccurate  as  the 
strength  of  the  electric  field  increases.  Unfortunately,  electric  field  strengths  generally  mcrease  as  devices  get 
smaller,  with  the  BSEs  becommg  essentially  worthless  for  predicting  the  performance  of  sub-micron  devices. 

We  have  recently  used  ringular  perturbatbn  techniques  to  analyse  the  Boltzmann  equatbns  for  semiconduc¬ 
tors.  This  analysis  shows  that  the  BSEs  are  an  asymptotically  correct  "outer  solution”  to  the  full  Boltzmann 
equations  for  moderate  electric  field  strengths.  However,  as  devices  become  smaller  and  electric  field  strengths 
increase,  the  Boltzmann  equatbns  enter  the  strong  force  regime,  which  renders  the  BSEs  incorrect.  In  this 
regime  our  analysis  yields  a  new  system  of  equations  in  which  the  drift  velocities  and  diffusivities  are  nonlin¬ 
ear  functions  of  the  electric  field.  Unlike  the  BSEs,  we  expect  these  new  High  Field  Semiconductor  Equations 
to  remain  accurate  regardless  of  the  electric  field  strength.  Remarkably,  these  equatbns  show  that  Einstem’s 
relation  is  satisfied  m  the  high  field  regime,  even  though  the  system  responds  nonlmearly  to  the  electric 
forces. 

We  will  also  derive  the  3-dimensbnal  anabgue  of  the  high  field  equatbns,  and  analyze  the  asymptotic  and 
quantum  mechanical  corrections  to  these  equations. 

3-3:30  pm  Leonid  V.  Kalachev  Some  applications  of  asymptotic  methods  in 

Moscow  State  University  semiconductor  device  modeling 

Abstract:  The  numerical  treatment  of  the  drift-diffusion  model  using  a  Gununel-type  iteration  process  often 
involves  the  soluton  of  nonlmear  equations  by  the  Newton  iteration  methods.  The  speed  of  convergence 
and  sometimes  the  convergence  itself  depend  crucially  on  the  right  choice  of  the  initial  iterate.  Asymptotic 
methods  show  their  effectiveness  when  used  to  construct  such  initial  iterates  as  asymptotic  approximations 
for  solutions  of  some  of  the  drift-diffusion  equations  when  a  small  parameter  enters  the  system.  We  discuss  the 
construction  of  the  asymptotics  for  the  singularly  perturbed  Poisson  equation  in  a  rectangular  semiconductor 
device.  Explicit  expressions  for  the  terms  of  the  asymptotic  approximation  significantly  provide  a  way  to 
use  asymptotics  m  numerical  processes  and  to  reduce  the  time  for  simulations. 

We  also  discuss  the  asymptotic  analysis  of  Gunn  effect  in  a  bounded  domain  modeling  the  Gunn  diode.  The 
asymptotic  approach  allows  one  to  construct  the  approximation  of  the  moving  pulse-type  solution  and  to 
obtain  some  characteristics  of  the  device. 
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We  give  a  short  survey  of  some  other  results  concerned  with  applications  of  asymptotic  methods  to  semicon* 
ductor  device  modeling. 

IMA  Fostdoc  Seminar 

4:00  pm  Chao*Nien  Chen  Bifurcation  for  prescribed  mean  curvature 

IMA  problems 

Ahsiraci:  We  discuss  bifurcation  and  existence  results  for  a  class  of  prescribed  mean  curvature  problen;s 
which  admit  zero  as  a  trivial  solution. 

The  SEMINAR  meets  in  Vincent  HaU  570 


Week  3  *•  Quantum  EifeCi..^  (\  week),  July  31  -  Aug.  2,  1991 
(Farouk  Odeh  and  Bill  Coughran  will  be  in  residence  and  in  charge  of  weeks  2  and  3) 

As  the  dimensions  of  modern  devices  decrease  to  ultrasmall  kvels,  quantum  effects  begin  to  play  a  major 
role.  Examples  are  the  quantisation  of  states  in  confined  regions  -  thin  channels,  quantum-wires  etc.  -  and 
tunneling  phenomena  in  resonant  diodes.  This  involves  a  coupling  of  electromagnetics  directly  to  quantum 
mechanics,  i.e.  to  eigenfunctions  of  the  Schroedmger  equation.  New  transport  models,  such  as  the  Wigner 
peeudo-diffeiential  equation  and/or  what  may  be  loosely  called  quantum-Bottsmann  have  to  be  applied  for 
an  adequate  description  of  ^quantum  devices’*.  Path  integrab,  pseudo-differential  operators  and  numerical 
techniques  for  eigenvalues  and  eigenfunctions  are  expected  to  provide  the  mathematical  tools  for  thb  area, 

Wednesday,  July  31 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 


9:30  am 

Harold  L.  Grubin 

Scientific  Research  Associates 

Small  devices  and  the  Boltzmann  equation  in 
quantum  mechanics 

10:30  am 

Coffee  Break 

Reception  Room  EE/CS  3-176 

llKKIam 

Christian  Ringhofer 

Arizona  State  University 

Numerical  methods  for  quantum  transport 
phenomena  in  semiconductor  devices 

2:00  pm 

H.C.  Liu 

National  Research  Council  Canada 

High  frequency  quantum  transport 

Abiir^ei:  A  scattering  model  is  formulated  for  quantum  transport  under  both  dc  and  ac  applied  vdtages. 
Analogous  to  the  scattering  model  for  the  time-independent  quantum  transport,  the  present  approach  views 
a  quantum  device  as  a  scattering  target  connected  to  electron  reservoirs  in  thermal  equilibrium.  Such  a 
model  is  therefore  applicable  only  to  the  extreme  quantum  cases,  e.g.,  tunneling  through  double-barriers 
and  ballistic  conduction  through  constrictions  or  quantum  point  contacts.  The  goal  is  to  investigate  (a)  how 
an  electron  interacts  with  a  high  frequency  ac  field  when  traversing  a  quantum  structure,  (b)  what  is  the 
frequency  scale  at  wnich  a  quantum  device  characteristics  deviate  from  its  low  frequency  values,  and  (c)  what 
is  the  device  high  frequency  characteristic,  e.g.,  device  impedance  as  a  function  of  frequency.,  A  physical 
insight  into  these  issues  U  of  great  relevance  to  the  frequency  limits  of  quantum  devices.  The  double-barrier 
resonanc  tunneling  Uode  and  the  bdlistic  constriction  are  discussed  as  examples. 

3-3:30  pm  Daniel  C.  Cole  The  thermodynamics  of  Casimir  forces  between 

IBM  GTD  conducting  parallel  plates 

Abstract  The  bulk  of  this  talk  is  related  to  a  recent  emerging  field  of  physics  called  cavity  quantum  elec¬ 
trodynamics.  However,  Casirnir  fcires,  as  well  as  the  underlying  van  der  Waals  forces,  also  play  a  role  in 
semiconductor  physics  and  devices,  which  Is  the  main  theme  of  this  conference.  This  role  is  briefly  discussed 
as  well  as  speculated  upon.  The  talk  then  turns  to  the  following  two  idealized  operations:  (1)  quasistatically 
displacing  two  parallel  conducting  plates  immersed  in  thermal  radiation,  and  (2)  slowly  changing  the  tem¬ 
perature  of  the  radiation.  By  imposing  the  condition  that  no  heat  must  flow  at  T  =  0  during  a  reversible 
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thermodynamic  operation,  the  functional  form  of  the  quantum  zero-point  radiaiton  spectrum  is  deduced. 
Ihming  to  thermodynamic  operations  on  the  plates  at  T  7^  0,  and  imposing  that  the  second  law  of  thermo¬ 
dynamics  must  hold,  then  results  in  a  generalized  derivation  of  Wien’s  displacement  law.  The  calculations 
involve  changes  in  singular  energy  expressions  and  cancellations  between  singular  forces. 

3:30  pm  Akerman  Hall  130D  Tour  of  Dan  Joseph’s  Lab 


We  will  show  experiments  on  water  lubricated  pipelining  with  bamboo  waves  and  corkscrew  waves,  two 
dimensional  cusped  interfaces  which  violate  L^lace’s  law,  rollers,  drafting,  kissing  and  tumbling  of  fluidized 
particles,  non  linear  stabilization  of  fingering  instabilities  in  porous  media,  drag  reduction  using  riblets.  The 
tour  probably  will  last  30  or  40  minutes. 


Thursday,  August  1 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 

9:30  am  Antonio  Gnudi  Some  i4>plication8  of  the  hydrodynamic 

University  of  Bologna  model  and  Boltzmann  transport  equation  in 

semiconductor  device  simulation 

Ab$iraci:  Phenomena  related  to  high  energy  carriers  in  semiconductors  play  an  important  role  in  the  design 
of  submicrometer  devices.  Some  models  have  been  devised  to  take  such  eflects  into  account  in  engineering 
simulation  tools  both  at  macroscopic  (Hydrodynamic  model)  and  microscopic  level,  by  means  of  a  more 
careful  examination  of  Boltzmann  Transport  Equation  (BTE). 

In  this  talk  we  review  some  applications  of  the  hydrodynamic  model  to  the  simulation  of  MOSFET’s,  with 
emphasis  on  the  ways  impact  ionisation  can  be  incorporated  in  the  model  in  a  physically  sound  yet  simple 
way.  Comparison  with  measured  substrate  currents  will  be  discussed. 

In  order  to  extract  nK>re  detailed  information  from  BTE,  a  method  will  be  outlined  based  on  expansion 
of  the  distribution  function  in  Spherical  Harmonics  in  momentum  space.  This  method  is  applied  to  both 
the  homogeneous  transport  problem,  with  comparison  to  Monte  Carlo,  and  to  a  one-dimensional  bipolar 
transistor. 

Joint  work  with  Farouk  Odeh. 

10:30  am  Coffee  Break  Pjeception  Room  EE/CS  3-176 

11:00  am  Erasmus  Langer  Numerical  simulation  of  MOS  transistors 

Technical  University  Vienna 

Abstract:  This  contribution  is  intended  to  review  the  international  state-of-the-art  in  numerical  simulation  of 
MOS  devices.;  Much  emphasis  is  laid  on  the  discussion  of  recent  refinements  to  carrier  transport  models,  e.g. 
drift-diffusion  model,  enhanced  drift-diffusion  equations,  hydrodynamic  model,  and  Monte  Carlo  simulation. 
Adequate  models  for  the  physical  parameters  are  reported  with  suitable  parameter  values,  e.g.  carrier  mobil¬ 
ities  taking  into  account  the  various  scattering  mechanisms,  and  carrier  generation-recombination  including 
impact  ionization.  Examples  are  presented  for  two  different  types  of  MOS  devices:  on  the  one  hand,  three- 
dimensional  simulation  results  of  a  miniaturized  MOS  transistor  are  discussed  which  have  been  obtained  by 
our  simulator  MINIMOS  5.0,  and  on  the  other  hand,  simulation  results  concerning  a  power  MOS  transistor 
are  shown  which  has  been  investigated  by  our  device  simulation  program  BAMBI  2.1.^ 

2:00  pm  B^^rnardo  Cockburu  Discontinuous  approximations  for  the  drift- 

U diversity  of  Minnesota  diffusion  and  the  hydrodynamic  models 

Abstract:  A  fully  parallelizable  method,  initially  devised  for  numerically  solving  nonlinear  conservation 
laws  in  several  ^pace  dimensions,  is  extended  to  the  drift-diffusion  and  to  the  hydrodynamic  models.  The 
method  is  con  bined  with  a  mixed  finite  eknient  that  provides  directly  an  approximation  of  the  electric 
field.  Stability  results  and  error  estimates  are  obtained  for  a  one-dimensional  version  of  the  drift-diffusion 
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model.  Preliminary  numerical  results  showing  the  performance  of  the  methods  are  shown  for  both  the 
(one-dimensional)  drift-diffusion  and  hydrodynamic  models. 

Joint  work  with  J.  Jerome  and  I.  IViandaf. 


EViday,  August  2 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 
This  day  b  set  aside  for  contributed  talks  and  dbcussions 

10:30  am  Coffee  Break  Reception  Room  EE/CS  3  ^76 


Week  4  -  Circuit  Analysb,  August  S-9,  1991 
(Jacob  White  will  be  in  residence  and  in  charge  of  week  4) 

The  high  cost  of  integrated  circuit  fabrication,  both  in  time  and  materials,  forces  demgners  to  simulate  their 
complete  designs  in  as  much  detail  as  possible.  Topically,  the  approach  of  choice  u  to  construct  a  differential 
algebraic  (DAE)  system  that  describes  the  circuit,  and  solve  that  system  numerically  for  a  ipven  set  of 
inputs.  More  specifically,  the  circuit  u  assumed  to  be  an  interconnection  of  lumped  elements  (e.g.  resutors, 
capacitor,  inductors  and  transbtors),  characterized  by  a  relation  between  the  terminal  voltages  uid  temunal 
currents.  The  equation  system  that  describes  such  a  network  b  then  constructed  using  the  Kirchoff  voltage 
and  current  laws,  and  this  leads  to  an  index  1  system  of  differential-algebraic  equations  where  the  state 
variables  are  not  always  readily  identified.  In  general,  the  system  of  equations  generated  from  most  circuits 
b  very  stiff,  even  if  the  algebraic  equations  are  not  present.  Abo.  the  system  u  quite  nonlinear,  and  thu 
makes  finding  an  initial  solution  which  satisfies  the  initial  conditions  and  the  algebraic  constraints  rather 
diflScult. 

The  circuit  simubtion  probbm  can  be  broken  into  two  classes:  solving  large  iuitial  value  problems,  with 
as  many  as  10,000  unlmowns,  efficwntly;  and  computing  periodic  or  quasi-periodic  solutions  for  sm^er 
problems,  with  10  to  100  unknowns.  For  solving  the  large  initial  v^ue  differential-algebraic  systems,  the 
topics  of  interest  are  techniques  for  stiff  systems;  tediniques  for  solving  large  nonlinear  systems  where  the 
Jacobian  can  be  much  more  than  one  rank  deficient;  multirate  integration  methods;  and  iterative  tech¬ 
niques  for  very  nonsymmetric  problems.  For  computing  periodic  or  quasi-periodic  solutions  to  these  DAE 
systems,  topics  of  current  interest  are  exbtence  and  uniqueness  of  solutions,  high  order  methods  for  very 
stiff  problems,  spectral  methods,  mixed  time-frequency  methods  and  other  types  of  coUocatiim  methods. 

Monday,  August  5 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 
9:30  am  Welcoming  Remarks  Conference  Hall  EE/CS  3-180 

9:40  am  Sally  Liu  Ir  Kishore  Singhal  Circuit  simulation  overview 

AT&T  Bell  Labs 

Abstract:  Recent  progress  in  circuit  simulation  will  be  reviewed,  bottlenecks  in  applying  current  state-of- 
the-art  circuit  simulation  capabilities  to  VLSI  and  high-performance  design  methodologies  will  be  presented. 
We  will  also  identify  some  open  issues  to  facilitate  discussions  at  the  workshop. 

10:40  am  Coffee  Break  Reception  Room  EE/CS  3-176 

11:00  am  Linda  R.  Petzold  On  the  numerical  solution  of  stiff  systems  of 

University  of  Minnesota  differential-algebraic  equations 

Abstract:.  In  recent  years  there  has  been  much  progress  in  the  analysis  and  development  of  numerical  methods 
and  software  for  the  solution  of  differential-algebraic  equation  (DAE)  systems.  However,  many  of  these 
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results  implicitly  assume  that  the  underlying  differential  system  is  nonstiff.  In  thb  lecture,  we  outline  recent 
developments  in  theory,  methods  and  software  which  have  implications  for  the  solution  of  stiff  DAEs.  In 
particular,  we  give  conditions  under  which  a  DAE  is  well-conditioned,  and  examine  formulations  of  the  DAE 
which  preserve  the  stability  and  which  lead  to  stable  discretizations.  For  initial  and  boundary  value  DAEs 
requiring  symmetric  discretizations,  we  describe  a  class  of  collocation  methods  which  overcome  stability 
and  accuracy  limitations  of  previously-defined  symmetric  methods.  Finally,  for  large-scale  stiff  initial  value 
DAEs,  we  discuss  the  status  of  available  software. 

2:00  pm  Bob  Melville  Globally  convergent  homotopy  methods  for  the 

AT&T  Bell  Labs  operating  point  and  steady-state  problems 

Ahsiraei:  Robust  computation  of  the  DC  operating  point(8)  of  an  integrated  circuit  is  crucial  for  simulation. 
An  operating  pomt  is  of  interest  in  its  own  right,  and  is  needed  as  a  starting  point  for  transient  analysis 
and  as  a  bias  point  for  small-signal  or  noise  analysis.  The  operating  point  problem  may  be  formulated  as 
a  system  of  non-linear  equations  to  be  solved  for  a  sero,  in  which  the  unknowns  are  node  voltages  and/or 
branch  currents.  In  some  cases,  the  system  of  equations  may  have  more  than  one  solution,  corresponding  to 
multiple  possible  bias  solutions  for  the  circuit  which  they  model. 

The  folbwing  three  methods  seem  to  be  the  most  popular  approaches  to  the  conq)utation  of  an  operating 
point  in  SPICE-class  programs: 

1.  Norm-reducing  variations  Newton’s  method; 

2.  Transient  methods  which  simulate  the  dynamics  of  the  circuit  as  the  power  supply  comes  on; 

3.  Numerical  continuation  which  tracks  the  state  of  the  circuit  for  different  values  of  the  power  supply 
voltage  (no<  as  a  dynamical  system  as  in  (2)  above). 

Our  presentation  makes  a  case  for  homotopy  (continuation)  parameters  other  than  supply  voltage.  In  par¬ 
ticular,  for  bipolar  networks,  the  use  of  transistor  current  gain  (forward  and  reverse  alpha)  as  a  homotopy 
parameter  gives  excellent  results.  Concepts  from  circuit  theory  are  used  to  establish  a  cotrciviiy  coniiiion  so 
that  the  sero  set  of  our  various  homotopiei  must  be  bounded.  However,  the  sero  set  can  still  exhibit  unde¬ 
sirable  singularities  such  as  bifurcations  and  cusps.  Such  singularities  can  be  eliminated  by  the  introduction 
of  an  appropriate  rsndomtjrsbon  vector.  Mathematical  justificatbn  for  this  procedure  will  be  given  based 
on  a  version  of  Sard’s  theorem. 

Timing  results  on  a  set  of  benchmarks  support  our  claim  that  genuine  global  convergence  can  be  eqjoyed  at 
a  computing  cost  which  is  easily  within  a  factor  of  five  of  other,  less  widely  convergent  techniques. 

The  steady-state  response  of  a  non-linear  network  driven  by  a  single  periodic  source  can  be  obtained  as  the 
solution  of  a  two-point  boundary  value  problem  on  a  discrete  mesh  of  time  points.  In  fact,  the  time-point 
equations  can  be  viewed  as  operating  point  equations  of  a  DC  network.  A  coercivity  result  will  be  given  for 
these  equations  which  lets  us  claim  global  convergence  of  a  homotopy  method  for  their  solution. 

Time  permitting,  we  will  discuss  multiple  solutions  of  operating  point  equations  and  qualitative  stability 
analysis  of  the  solutions. 

Various  parts  of  the  work  to  be  presented  have  been  joint  with  Ljiljana  IVajkovic  at  BellCore,  and  Mike 
Fang,  Peter  Doyle,  and  of  Jeff  Lagarias  at  AT&T. 

4:00  pm  Vincent  Hall  502  IMA  Tea  (and  more!) 

(The  IMA  Lounge) 


9.30  am 


Tuesday,  August  6 

Unless  otherwise  stated,  the  talks  today  cure  in  Conference  Hall  EE/CS  3-180 

Serial  and  parallel  waveform  relaxation 
algorithms 


Albert  Ruehli 
IBM 


10:30  am  Coffee  Break 


Reception  Room  EE/CS  3-176 
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11:00  am  Donald  J.  Rose  Newton-iterative  and  iterative-Newton  methods 

Duke  University  for  nonlinear  systems 

Abstract:  We  discuss  two  generic  methods  for  solving  nonlinear  algebraic  or  operator  systems  of  equations. 
Newton-iterative  methods  solve  successive  linearized  systems  by  an  inner  iterative  method,  i.e.,  Newton- 
Gauss^Seidel.  Iterative-Newton  methods,  as  in  Guass-Seidel-Newton,  solve  smaller  sets  of  nonlinear  equations 
approximately,  usually  using  one  step  of  Newton’s  methods.  We  compare  these  methods  with  an  eye  toward 
computational  efficiency;  in  particular,  memory  requirements,  function  evaluations,  arithmetic  operations, 
and  ease  of  implementation.  Circuit  and  device  simulation  motivate  our  discussion. 

2:00  pm  A.  Lumsdaine  Gradient-descent  acceleration  for  waveform 

MIT  relaxation 

Abstract:  A  new  algorithm  for  accelerating  dynamic  iterations  for  DAE  systems  is  developed.  The  algorithm 
is  based  on  extend*  the  generalized  coiyugate-residua!  algorithm  for  matrix  problems  to  function  spaces. 
Experimental  results  for  the  method  used  to  solve  time-dependent  drift-diffusion  plus  Poisson  equation  is 
given  to  demonstrate  the  effectiveness  of  the  acceleration. 

3:00  pm  Ben  Leimkuhler  Estimating  waveform  relaxation  convergence 

University  of  Kansas 

Abstract:  The  waveform  relaxaticHi  method  is  considered  as  a  method  for  decoupling  systems  of  ODEs  or 
DAEs  for  parallel  implementation.  Ultimately,  waveform  relaxation  converges  superlinearly  in  finite  intervals. 
In  the  early  sweeps,  this  is  typically  not  observed,  and  it  has  been  proposed  that  the  practical  convergence 
is  linear  until  the  late  stages  of  the  iteration.  Taking  this  viewpoint,  the  author  obtains  an  estimate  for 
a  window  length  wherein  waveform  relaxation  converges  like  w”  by  looking  in  the  Laplace  domain.  The 
results  are  stated  in  terms  of  a  quantity  called  the  speed  of  a  splitting  which  measures  the  extent  to  which 
subsystems  are  decoupled.  A  numerical  experiment  is  presented. 


Wednesday,  August  7 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3*180 

9:30  am  Larry  Pillage  Moment-matching  approximations  for 

University  of  Texas  linear(ized)  circuit  analysis 

Abstract:  Moment-matching  approximations  appear  to  be  a  promising  approach  for  linear  circuit  analysis 
in  several  application  areas.  Asymptotic  Waveform  Evaluation  (AWE)  uses  moment-matching  to  approx¬ 
imate  the  time-  or  frequency-domain  circuit  response  in  terms  of  a  reduced-order  model.  AWE  has  been 
demonstrated  as  an  efficient  means  for  solving  large,  stiff,  linear  circuits.  However,  since  it  is  based  upon 
moment-matching,  which  has  been  shown  to  be  equivalent  to  a  Fade’  approximation,  AWE  is  prone  to 
yielding  unstable  waveform  approximations  for  stable  circuits.  In  addition,  it  is  difficult  to  quantify  the 
time  domain  error  for  moment-matching  approximations..  We  address  the  issues  of  stability  and  accuracy  of 
moment-matching  approximations  as  they  apply  to  linear  circuit  analysis. 

10:30  am  Coffee  Break  Reception  Room  EE/CS  3-176 

11:00  am  Jacob  Wliite  Multipole- accelerated  3-D  capacitance  extraction 

MIT 

Abstract:  In  this  talk  we  will  describe  a  fast  method  for  computing  the  capacitance  matrix  of  a  general 
3-dimensional  geometry  of  ideal  conductors  in  free  space.  The  method  is  an  acceleration  of  the  standard 
boundary-element  approach  to  solving  the  first-kind  integral  equation  for  cohductor  surface  charge.  In  tins 
accelerated  approach,  the  dense  matrix  problem  generated  from  the  boundary-element  discretization  is  solved 
with  a  preconditioned  generalized  conjugate-residual  method,  where  a  multipole  approximation  is  used  to 
compute  the  iterates,  and  the  preconditioner  is  matched  to  the  multipole  algorithm  A  novel  adaptive  scheme 
for  the  3-dimensional  inultipole  algorithm  will  be  presented,  along  with  results  indicating  the  effectiveness  of 
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the  c^proach  on  engineering  problems.  Finally,  numerical  evidence  is  presented  indicating  that  the  number 
of  GCR  iterations  does  not  increase  with  increasing  geometric  complexity,  demonstrating  that  in  practice 
the  algorithm  is  order(iV),  where  iV  is  the  number  of  boundary  elements. 

2:00  pm  Olgierd  A*  Palusinski  Spectral  technique  in  simulation  of  MOS  circuits 

University  of  Arizona,  Ihcson 

Abstract:  Waveform  Relaxation  (WR)  is  an  attractive  way  to  reduce  CPU  time  needed  for  simulation  of 
transients  m  integrated  circuits.  One  of  the  problems  associated  with  WR  is  a  need  for  storage  of  intermediate 
solutions  for  coupling  variables  and  exchange  of  those  variables  between  sq>propriate  subcircuits.  Spectral 
technique  based  on  Chebyshev  polynomials  offers  the  most  compact  representation  of  variables  and  has 
many  features  facilitating  computation  and  exchange  of  data  which  is  necessary  in  WR.  The  technique  will 
be  briefly  described  and  some  results  of  application  to  MOS  circuits  ipven. 


Thursday,  August  8 

Unless  otherwise  stated,  the  talks  today  are  in  Conference  Hall  EE/CS  3-180 

9:30  am  Jiri  Vlach  Switched  networks 

University  of  Waterloo 

Abstract:  A  special  integration  method,  based  on  numerical  Laplace  inversion,  was  developed  for  the  analysis 
of  switched  networks.  The  method  can 

(a)  find  correct  time  domain  response  in  the  presence  of  Dirac  impulses 

(b)  find  consistent  initial  conditions  after  switching,  even  if  at  the  time  of  switching  the  conditions  are  inconsutent 
and  Dirac  impulses  may  be  present. 

(c)  Discover  not  only  the  absence  or  presence  of  the  Dirac  impulse,  but  also  its  area  and  perform  correct  switching 
of  other  elements  in  the  network. 

The  method  was  used  for  the  development  of: 

(a)  Program  for  analysis  of  periodically  switched  analog  networks  (with  all  linear  elements). 

(b)  Program  for  analysis  of  networks  with  internally  controlled  switches,  as  needed  for  analysis  of  switched  power 
supplies. 

The  lecture  will  explain  basis  of  the  integration  method  and  will  show  examples  solved  by  our  programs. 
10:30  am  Coffee  Break  Reception  Room  EE/CS  3-176 

11:00  am  Michael  £.  Henderson  Computing  periodic  solutions  and  their 

'  IBM,  Watson  Research  Center  parameter  dependence 

Abstract:  It  has  become  fairly  standard  to  compute  periodic  solutions  in  systems  of  ODE’s,  and  to  track  the 
solutions  and  their  bifurcations  as  some  parameter  in  the  equations  is  varied.  I  will  describe  the  equations 
used  and  some  of  the  considerations  which  arise  when  the  method  is  implemented.  Finally,  I  will  describe 
one  application,  where  these  techniques  have  been  used  to  compute  periodic  solutions  in  a  pair  of  coupled 
Josephson  junctions. 

2:00  pm  Ken  Kundert  Mixed  frequency-time  methods  for  multiperiodic 

Cadence  Design  Systems  systems 

Abstract:  Analog  circuits  like  switching  filters  are  typically  driven  with  a  clock  which  is  at  a  much  higher 
frequency  than  the  signals  the  circuit  is  intended  to  process.  Designers  of  such  circuits  are  often  interested 
in  steady-state  distortion  due  to  both  static  effects,  such  as  nonlinearities  in  the  capacitors,  and  dynamic 
effects,  such  as  the  charge  injection  during  MOS  transistor  switching  or  slow  operational  amplifier  settling. 
Steady-state  distortion  can  be  computed  using  the  circuit  simulation  program  SPICE,  but  this  approach  is 
computationally  very  expensive.  Specialized  programs  for  switched  capacitor  filters  can  be  used  to  rapidly 
compute  steady-state  distortion,  but  do  not  consider  dynamic  effects  In  this  talk  we  present  a  new  mixed 
frequency-time  approach  for  computing  both  steady-state  distortion.  The  method  is  both  computationally 
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efficient  and  includes  both  static  and  dynamic  distortion  sources.  The  method  has  been  implemented  in  a 
C  program,  NtUvrii,  and  results  from  several  examples  will  be  presented. 

3:00  pm  Werner  Liniger  On  sparse  factorisation  methods  for  Poisson's 

IBM  Yorktown  Heights  equation 

Abstmct  We  discuss  fast  algorithms  for  solving  Poisson’s  (or  other  elliptic)  equations  on  general  regions 
with  arbitrary  mixed  boundary  conditions.  These  algorithms  create  preconditioners  which  are  second-order, 
both  in  the  interior  of  the  region  and  near  the  boundary.  The  associated  iterative  scheme  has  interesting 
and  novel  convergence  properties. 

5:30  pm  S02  te  570  Vincent  Hall  Gismt  Hoagie  Party  (sandwiches,  salad, 

fruit,  desert) 

fViday,  August  9 

The  talk  today  is  in  Vincent  Hall  570 
9:30  am  Paul  Cox  Circuit  simulation  algorithnw 

Texas  Instruments 

Ab$trmet:  Circuit  simulators  are  one  of  the  most  widely  used  and  essential  tools  used  in  the  design  of 
integrated  circuits,  lb  satisfy  these  requirements  two  conflicting  requirements  nuist  be  sidisfled  'Hie  circuit 
sinnulation  must  provide  the  accuracy  and  robustness  required  for  critical  memory  and  anab*  .ireuit  designs 
and  it  must  be  efficient  enough  to  simulate  large  VLSI  circuit  designs.  We  have  been  involved  for  some  time 
in  the  development  and  refinement  of  the  circuit  simulator  SUPPLE  at  Texas  Instruimnts.  The  objective  of 
this  work  was  to  provide  a  circuit  simulator  with  all  of  the  capabilities  and  robustness  of  SPICE  but  with 
enhanced  performance  from  the  exploitation  of  nmilti-rate  behavior  and  parallel  processing. 

The  algorithms  developed  for  use  in  SUPPLE  have  improved  the  efficiency  to  simulation,  however,  a  number 
of  issues  need  to  be  resolved  in  order  to  support  the  design  of  larger  VLSI  designs.  These  issues  include: 
DC  convergence.  A  DC  solution  of  the  circuit  matrix  is  the  first  step  in  any  type  of  circuit  simulation. 
However,  no  efficient  noethods  exists  which  can  guarantee  a  DC  solution. 

Matrix  partitioning.  Efficient  methods  of  partitioning  need  to  be  developed  for  exploitation  of  multi-rate 
behavior  and  parallel  processing. 

IVansient  time  step  control.  New  algorithms  need  to  be  developed  to  provide  better  a  priori  truncation  error 
estimation. 

Matrix  solution  algorithms.  New  matrix  solution  techniques  are  required  which  can  provide  the  robustness 
of  direct  solution  methods  but  are  efficient  enough  to  be  used  on  very  large  sparse  matrix  systems. 

10:30  am  Coffee  Break  Vincent  Hall  502 


The  rest  of  this  day  is  set  aside  for  discussions 


t  CURRENT  IMA  PARTICIPANTS  "1 


NAME 


POSTDOCTORAL  MEMBERS  FOR  1990-91  PROGRAM  YEAR 

PREVIOUS/PRESENT  INSTITUTION 


Affouf,  Mahmoud 
Chen,  ChaoNien 
Choe,  Hi  Jun 
Fan,  Haitao 
Fehribach,  Joseph 
Filippas,  Stathis 
Firoozye,  Nikan 


New  Jersey  institute  of  Technology 

Indiana  University 

University  of  Kentucky 

Brown  University 

University  of  Alabama,  Huntsville 

Courant  Institute 

Courant  Institute 
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Mou,  Libin  H.  Rice  University 

Smereka,  Peter  Courant  Institute 


POSTDOCTORAL  MEMBERSHIPS  IN  INDUSTRIAL  MATHEMATICS  FOR  1990-91  YEAR 
NAME  PREVIOUS/PRESENT  INSTITUTION 


Dobson,  David  C. 
F^emi,  Emad 
Hoffend,  Thomas  R.  Jr.. 
Morokoff,  William 
Xu,  Yongshi 


Rice  University 

University  of  California,  Los  Angeles 
State  University  of  New  York,  Buffalo 
Courant  Institute 
University  of  Delaware 


LONG-TERM  VISITORS  IN  RESIDENCE 
4  Weeks  or  Longer 


DonatOi  Jerry 

Jun  l(9l)-Jun  30(92) 

Dreiin,  Yuri  A* 

IMA 

June  1  -  Nov  30 

FViedman,  Avner 

IMA 

Guo,  Jong-Sheng 

National  Tsing  Hua  University 

Jul  15~Aug  15 

Hu,  Bei 

University  of  Notre  Dame 

June  10  -  Aug  15 

Li,  Ibng 

Princeton  University 

Feb  10  -  July  20 

Littman,  Walter 

University  of  Minnesota 

Sept  1  -  June  30 

MiUer,  WiUaid,  Jr. 

IMA 

Mukherjee,  S. 

Mahatma  Gandhi  Degree  Coll. 

Aug  27  (90)  -  Aug  26,  (91) 

Ssmigieltki,  Jacek 

University  of  Virginia 

Jun  17  -  Aug  15 

TVan,  Victor  Ngoc 

UC  Irvine 

Jun  27  •  Aug  28 

von  Schwerin,  Reinholl 

Universitat  Augsburg 

May  8  -  Oct  31  (91) 

Wang,  Lihe 

Princeton  University 

Feb  10  -  July  20 

Weinberger,  Hans 

University  of  Minnesota 

SUMMER  PROGRAM  PARTICIPANTS  IN  RESIDENCE 

Aarden,  J. 

University  of  Ngmegen 

Jul  20-Aug  3 

Baccarani,  Giorgo 

University  of  Bologna 

Jul  22- Aug  2 

Bennett,  Herbert 

NIST 

Jul  21-Aug  2 

Biswas,  Rana 

Iowa  State  University 

Jul  16  -  Jul  17 

Blakey,  Peter 

Motorola  Corporation 

Jul  21-Aug  2 

Ebrucki,  Leonard 

Motorola  Corporation 

Jul  14-Jul  19 

Buergler,  Josef 

ETH  Zurich 

Jul  21  -  Aug  2 

Casey,  Michael 

University  of  Pittsburgh 

Jul  15-Aug  9 

Cercignani,  Carlo 

Politecnico  di  Milano 

Jul  21-Aug  3 

Cole,  Dan 

IBM  GPD 

Jul  24- Aug  1 

Cole,  Julian 

Rensselaer  Polytechnic  Institute 

Jul  14-Aug  2 

CougKran,  Jr.,  William 

AT&T  Bell  Ubs 

Jul  28-Jul  31 

Cox,  Paul 

Texas  Instruments 

Aug  3-Aug  9 

Degond,  Pierre 

Ecole  Polytechnique 

Jul  21-Aug  2 

Gaal,  Steven 

University  of  Minnesota 

Jul  15  -  Aug  9 

Gardner,  Carl 

Duke  University 

Jul  22-23,  29-31 

Gartland,  Chuck 

Kent  State  University 

Jul  2i-Jul  31 

Gerber,  Etean 

IBM 

Jul  14-Jul  19 

Giles,  Martin 

University  of  Michigan 

Jul  14-Jul  16 

Glodjo,  Arman 

University  of  Manitoba 

Jul  14-Aug  9 

Gnudi,  Antonio 

Universita  Degli  Studi  Di  Bologna 

Jul  30- Aug  2 

Grubin,^  Harold 

Scientific  Research  Associates 

Jul  30- Aug  2 

Hagan,  Patrick 

Los  Alamos  National  Lab 

Jul  2 1-24 .,28-31 
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Hamaguchi,  Satoshi 
Henderson,  Mike 
Jerome,  Joseph  W. 
Johnson,  Michael 
Kalachev,  Leonid 
Kerkhoven,  Thomas 
King,  John 
Kundert,  Ken 
Langer,  Erasmus 
Law,  Mark 
Leimkuhler,  Ben 
Liniger,  W. 

Utt,  H.C. 

Uu,  Sally 
Liu,  Xu'Dong 
Lloyd,  Peter 
Loj^,  Robert 
Lumsdaine,  Andrew 
Makohon,  Richard 
Meinershagen,  Bemdt 
Melville,  Robert 
O'Malley,  Robert  E. 
Odeh,  Fuouk 
Palusinski,  O. 

Perline,  Ron 
Petsold,  Linda  R. 
Pidatdla,  Rosa  Maria 
Pillage,  Larry 
Please,  Ciolin 
Poupaud,  Frederic 
Royna,  Luis 
Rkhardaon,  Walter 
Ringhofer,  Christian 
Rose,  Donald  J. 
Rudan,  Massimo 
Ruehli,  Albert 
Schmeiser,  Christian 
Seidman,  Tom 
Sever,  Michael 
Singhal,  K. 

So,  Wasin 
Souissi,  Kamel 
Strojwas,  Andre 
Suto,  Ken 
Szmolyan,  Peter 
Tang,  Henry 
Thomann,  Enrique 
Venturino,  Ezio 
Vlacn,  Jeri 
Ward,  Michael 
White,  Jacob 
Wrzosek,  Darek 
Young,  Richard  A 


IBM 

IBM 

Northwestern  University 
IBM 

Moscow  State  Univernty 
University  of  Illinou,  Urbana 
University  of  Nottingham 
Cadence  Design  Systems 
Technical  U.  Vienna 


Jul  14-Jul  19 
Aug  4-Aug  9 
Jul  21-Jul  26 
Jul  14-Jul  19 
Jul  1&-Aug  10 
Jul  23-Aug  1 
Jul  14-Jul  25 
Aug  4-Aug  9 
Jul  21-Aug  2 


University  of  Florida  Jul  14-Jul  19 

University  of  Kansas  Jul  14-Aug  9 

IBM  Aug  3-Aug  9 

National  Research  Council,  Ottawa  Jul  25-Aug  3 
AT&T  BeU  Labs  Aug  4-Aug  9 

UCLA  Jul  14-Jul  21 

AT&T  BeU  Ubs  Jul  14-Jul  IS 

Motorola  Jul  14-Jul  23 


Mrr  Aug  3-Aug  9 

University  of  Portland  Jul  14-Aug  2 

Technischen  Hodisdiule  Aadien  Jul  21-Aug  2 

AT&T  BeU  Labs  Aug  3-Aug  9 

Rensselaer  Polytechnic  Institute  Jul  21-Aug  3 

IBM  Jul  21-Jul  31 


University  of  Arizona 
Drexel  University 
University  of  MinnesrAa 
CStta*  Univerntaria,  Italy 
University  of  Texas 
Southhampton  University 
University  of  Nice 
IBM 


Aug  3-Aug  9 
Jul  14-Aug  9 
Jul  24-Aug  9 
Jul  13-Aug  10 
Aug  3-Aug  9 
Jul  16-Jul  27 
Jul  21-Jul  31 
Jul  23-Jul  28 


University  of  Texas  at  San  Antonio  Jul  14-Aug  9 


Arizona  State  University 
Duke  University 
University  of  Bologna 
IBM 

TU-Wien-Austria 

U.  of  Maryland-Baltimore  County 

Hebrew  University 

AT&T  BeU  Ubs 

IMA 

IBM 

Carnegie  Mellon  University 
Tohoku  University 
TU-Wien-Austria 
IBM 


Jul  28-Aug  1 
Aug  4-Aug  9 
Jul  20-Jul  30 
Aug  4-Aug  9 
Jul  21-Aug  2 
Jul  21-Aug  2 
Jul  21-Aug  3 
Aug  4-Aug  9 
Jul  16  -  Aug  31 
Jul  27- Aug  3 
Jul  14-Jul  19 
Jul  14-Jul  19 
Jul  21-Aug  2 
Jul  21-Aug  2 


Oregon  State  University 
University  of  Iowa 
University  of  Waterloo 
Stanford  University 


Jul  20-Jul  29 
Jul  14-Aug  9 
Aug  3-Aug  9 
Jul  25-Jul  29 


MIT 

University  of  Warsaw 
University  of  Portland 


Aug  3-Aug  9 
Jul  14-Jul  31 
Jul  14-Aug  2 
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